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Abstract 

The  Elnited  States  Air  Eoree  is  partnering  with  the  Einited  States  Army  as  well  as 
allied  nations  to  develop  a  revolutionary  advanee  in  logistical  support  known  as  the  Joint 
Precision  Air  Drop  System  (JPADS).  The  focus  of  this  study  is  to  develop  a  process  to 
quantitatively  analyze  system  sensitivities  to  various  types  of  weather  inputs  and  the 
corresponding  effect  on  system  accuracy.  Weather  balloons  were  used  to  provide 
representative  “truth”  to  which  forecast  weather  could  be  compared.  Each  data  type  was 
fed  into  the  JPADS  Mission  Planner  to  produce  navigation  points  which  could  then  be 
compared  statistically.  The  process  was  tested  on  a  limited  data  set  to  provide  a  first  look 
at  the  variables  of  forecast  resolution  and  “lead-time.”  Initial  results  indicate  best  system 
accuracy  is  achieved  for  lowest  forecast  resolution  (i.e.,  45  km  vs.  5  km  data)  and 
shortest  lead-time  (i.e.,  <12  hrs  vs.  >12  hrs).  This  result  will  not  only  allow  for  better 
accuracy  of  JPADS,  but  also  reduce  bandwidth  and  transmission  time  necessary  to  send 
weather  forecast  data  to  the  warfighter. 
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A  SYSTEMS  ENGINEERING  APPROACH  TO  ANALYZING  WEATHER  INPUT 


SENSITIVITIES  OE  THE  JOINT  PRECISION  AIR  DROP  SYSTEM 


I.  Introduction 


Background 

In  the  modern  world  of  precision  engagement  with  weapons  such  as  the  Joint 
Direct  Attack  Munition  (JDAM)  allowing  for  accuracy  measured  in  feet,  the  airdrop 
community  has  had  to  soldier  on  with  low  precision  tactics  and  techniques  that  would  be 
recognizable  to  their  Vietnam  era  counterparts.  The  Joint  Precision  Air  Drop  System 
(JPADS)  is  a  family-of-systems  developed  to  bring  the  sort  of  precision  capability  found 
in  Global  Positioning  System  (GPS)  guided  munitions  to  the  airlift  community.  As  such, 
JPADS  is  often  touted  as  “the  JDAM  for  logistics.” 

The  purpose  of  this  thesis  is  to  investigate  the  impact  of  weather  data  inputs  on 
the  accuracy  of  JPADS,  specifically  the  JPADS  Mission  Planner  (JPADS-MP)  and  the 
navigation  outputs  it  creates  from  this  weather  data.  Despite  the  conceptual  similarities 
between  JDAM  and  JPADS,  the  airdrop  mission  poses  problems  for  accuracy  not  faced 
in  the  precision  guided  munition  mission.  A  JDAM  class  munition  falls  at  high  velocity 
through  the  atmosphere  and,  while  guided,  still  follows  a  relatively  ballistic  trajectory. 

As  such,  the  precision  of  JDAM  type  munitions  is  not  greatly  affected  by  the  weather 
they  pass  through  between  the  launching  aircraft  and  their  target.  The  JPADS  chute 
loads  are  quite  different  as  they  truly  do  fly  as  a  paraglider.  Since  they  are  unpowered, 
proper  energy  management  during  their  decent  is  critical  in  hitting  their  designated  Point- 
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of-Impact  (PI).  In  order  to  achieve  the  desired  level  of  precision,  a  Guided  Parachute 
requires  knowledge  of  the  state  of  the  atmosphere  in  which  it  will  fly.  This  thesis 
examines  the  current  methods  used  for  weather  data  ingestion  by  JPADS  and  determines 
best  practices  for  the  as-is  system.  It  will  conclude  with  recommendations  for  further 
research  to  develop  an  improved  to-be  system. 

Scope 

Even  as  this  thesis  is  being  written,  JPADS  is  evolving.  At  present,  it  is 
undergoing  the  Joint  Military  User  Assessment  stage  of  its  Advanced  Concept 
Technology  Demonstration  (ACTD)  testing  program  at  the  US  Army’s  Yuma  Proving 
Ground  (YPG).  Despite  the  continued  evolution,  JPADS  is  already  operating  in  the 
combat  theatre.  The  fact  that  JPADS  is  already  in  use  serves  to  focus  the  domain  of  this 
thesis.  This  research  is  the  product  of  techniques  from  weather  forecasting,  operations 
research,  and  systems  engineering.  Any  one  of  these  fields  can  find  a  rich  source  of 
problems  for  study  in  the  JPADS  program.  However,  there  is  a  pressing,  operational 
question  at  hand.  As  the  system  stands  today,  in  the  field,  how  do  we  make  it  as  accurate 
and  precise  as  possible?  Other  organizations  are  already  pursuing  studies  of  the 
Guidance,  Navigation,  and  Control  algorithms  for  JPADS.  That  leaves  the  question  of 
weather  impacts. 

Problem  Statement 

Airdrop  operators  require  an  evaluation  of  the  sensitivity  of  the  JPADS-MP  to 
weather  inputs.  To  quantify  this  sensitivity,  it  is  necessary  to  first  identity  what  weather 
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products  are  used  by  JPADS-MP  and  what  they  are  used  for.  To  aehieve  a  manageable 
scope,  this  study  will  focus  on  two  weather  produets  usable  by  the  JPADS-MP;  Air  Foree 
Weather  Agency  (AFWA)  Forecasts  and  Weather  Balloons.  The  immediate  goal  is  to 
statistieally  eompare  the  various  types  of  foreeasts  generated  by  AFWA  to  aetual  weather 
sampled  by  weather  balloons  and  thus  determine  the  best  operational  practice.  Of  even 
greater  value,  though,  is  the  proeess  that  will  be  developed  to  aehieve  this  goal,  as  it  will 
eontinue  to  be  useful  to  the  program  as  an  analytieal  proeess  beyond  this  initial  research. 

Research  Objectives,  Questions,  and  Hypotheses 

Research  Objective. 

The  objeetive  of  this  study  is  to  analyze  the  weather  sensitivities  of  the  Joint 
Preeision  Air  Drop  System.  To  do  this  requires  the  development  of  a  standardized, 
statistieally  sound,  method  of  eomparing  weather  inputs  to  the  JPADS-MP.  This  researeh 
will  then  use  this  proeess  to  perform  an  initial  analysis  to  answer  the  Researeh  Question. 

Research  Question. 

This  researeh  will  perform  an  initial  analysis  of  AFWA  weather  foreeasts  to 
determine  whieh,  if  any,  provide  better  planning  aeeuraey  for  the  airdrop  mission  when 
used  as  input  to  the  JPADS-MP. 

Investigative  Questions. 

1 .  How  does  weather  affeet  JPADS?  (i.e.,  how  does  the  JPADS-MP  ingest  and 
use  weather  data?  What  are  the  outputs?) 

2.  How  does  the  JPADS-MP  use  weather  data  to  generate  navigation  outputs? 
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3.  How  can  the  navigation  outputs  from  the  JPADS-MP  be  eonverted  to  a 
statistieally  eomparable  format? 

4.  What  different  types  of  AFWA  forecasts  are  available? 

5.  How  are  weather  balloons  made  available  to  the  JPADS-MP? 

6.  What  statistieal  tests  and  tools  can  be  used  to  analyze  the  weather  sensitivities 
of  JPADS? 

Hypotheses. 

A  key  objeetive  of  this  study  is  to  apply  statistics  in  order  to  get  a  quantitative 
understanding  of  how  JPADS  is  sensitive  to  weather  inputs.  As  previously  mentioned, 
JPADS  is  already  operating,  both  in  test  and  in  the  field,  and  is  doing  so  with  eertain 
qualitative  assumptions  about  the  best  practiees  to  use  regarding  weather  inputs.  JPADS 
has  already  made  great  strides  in  aeeuraey,  but  more  is  desired.  To  get  there,  qualitative 
assumptions  must  give  way  to  quantitative  results. 

The  statistieal  tests  used  in  this  study  are  relatively  simple.  In  general,  the  Null 
Hypothesis  (Ho)  will  be  that  the  given  distribution  eannot  be  rejeeted;  while  the  Alternate 
Hypothesis  (Ha)  will  be  that  the  given  distribution  is  rejeeted.  All  statistieal  tests  in  this 
study  will  be  performed  at  a  =  0.05. 

Methodology  in  Brief 

A  standardized  mission  will  be  used  in  the  JPADS-MP,  with  the  different  weather 
produets  being  used  to  generate  navigation  outputs.  These  outputs  will  then  be  eonverted 
into  a  eommon  Northing  vs.  Easting  error  format  whieh  is  funetionally  like  the  Miss 
Distanee  eharts  used  by  the  JPADS  Guided  Paraehute  systems.  The  data  will  then  be 
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subjected  to  statistical  examination  to  determine  Goodness-of-Fit  for  Bivariate  Normality 
(a  typical  distribution  for  this  type  of  data).  Finally,  the  means  and  variances  of  the 
different  data  groups  were  compared  to  identify  the  best  weather  forecast  type  for  use. 

Document  Structure 

Chapter  2  is  a  literature  review  which  will  provide  a  more  in  depth  discussion  of 
the  topics  introduced  in  Chapter  1 .  The  chapter  will  begin  with  an  abbreviated  history  of 
airdrop  then  progress  to  a  brief  description  of  the  JPADS  system,  as  well  as  provide  some 
insight  into  weather  forecasting  and  weather  products  used  by  JPADS-MP.  Chapter  3 
details  the  means  by  which  the  research  was  accomplished.  It  includes  further 
information  on  how  the  weather  data  was  converted  into  a  useful  format  for  statistical 
analysis  as  well  as  the  details  of  that  analysis.  Chapter  4  reports  the  results  of  the 
analysis  and  Chapter  5  provides  a  summary  of  the  overall  research  effort  as  well  as 
presents  avenues  for  further  research. 

Limitations 

The  research  documented  in  this  thesis  is  limited  to  the  analysis  of  historical  data 
rather  than  a  fully  designed  original  experiment.  As  such,  the  analysis  in  this  research 
must  use  data  which  was  intended  for  other  purposes.  In  this  case,  not  all  potentially 
available  forecasts  were  recorded.  While  there  are  sufficient  data  points  to  extract 
statistical  significance,  care  must  be  taken  in  interpreting  the  results  so  as  not  to  over 
generalize  beyond  what  the  data  supports. 
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The  weather  balloons  that  were  used  as  a  basis  of  comparison  for  the  forecast  data 
were,  of  course,  actually  launched  to  support  the  aircraft  operations  of  the  JPADS  ACTD. 
This  makes  it  necessary  to  check  the  data  for  unanticipated  correlations.  Additionally,  it 
is  worth  noting  that  the  primary  research  question,  i.e.  the  accuracy  of  various  forecast 
products,  was  originally  raised  by  AFWA  Det  3.  It  was  decided  to  use  the  JPADS-MP  as 
an  analysis  tool  rather  than  analyze  the  various  weather  data  products  directly.  This 
indirect  method  was  chosen  for  two  reasons.  First,  the  JPADS-MP  must  perform  an 
internal  analysis  in  order  to  generate  navigational  outputs;  and  second  by  using  the  actual 
mission  planning  tool  the  operators  use,  the  warfighter  is  assured  of  a  result  with 
immediate  operational  application. 

Finally,  all  forecast  data  used  for  this  thesis  was  collected  for  YPG  and  used  the 
Penn  State  University/National  Center  for  Atmospheric  Research  Mesoscale  Model  5 
(MM5)  forecast  model.  This  analysis  will  need  to  be  reaccomplished  when  AFWA 
changes  from  MM5  to  the  Weather  Research  and  Forecasting  (WRF)  model.  Time 
limitations  prevented  attempting  to  gather  data  from  areas  other  than  YPG.  It  is  therefore 
worthwhile  to  use  caution  in  applying  the  results  of  this  research  to  other  locales  before 
additional  data  can  be  reviewed.  The  methodology  used  in  this  thesis  will  allow  for  such 
additional  analysis  with  ease.  This  is  an  advantage  of  using  the  JPADS-MP  as  an 
analysis  tool. 
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II,  Literature  Review 


Overview 

The  Joint  Precision  Air  Drop  System  is  intended  to  address  several  recognized 
capability  gaps.  It  is  a  family  of  systems  that  includes,  but  is  not  limited  to,  the  JPADS 
Fly-Away  Kit  and  several  candidate  guided  parachute  systems.  This  chapter  will  begin 
with  a  brief  discussion  of  the  historical  environment  that  led  to  JPADS.  It  will  then 
progress  to  a  description  of  the  systems  that  comprise  JPADS.  This  will  be  cursory  as 
JPADS  is  well  covered  in  other  documents  and  is  not  the  actual  focus  of  this  thesis. 
Attention  will  be  given  to  aspects  of  the  JPADS-MP  which  were  of  specific  use  in  this 
thesis.  The  chapter  will  conclude  with  a  review  of  the  weather  data  types  and  formats 
used  in  this  research. 

Historical  Background 

On  16  January,  1784,  an  American  living  abroad  in  France  penned  a  letter  to  a 
friend  concerning  a  revolutionary  technology  he  had  recently  observed.  The  technology 
in  question  was  a  balloon  capable  of  lifting  two  men  into  the  air.  This  American  saw 
more  than  a  mere  curiosity  in  the  balloon.  In  fact,  he  had  an  extraordinarily  prescient 
vision  of  what  would  stem  from  the  invention. 

On  that  day  he  wrote: 

. .  .where  is  the  prince  who  can  afford  so  to  cover  his  country  with  troops 
for  its  defense,  as  that  ten  thousand  men  descending  from  the  clouds  might 
not  in  many  places  do  an  infinite  deal  of  mischief,  before  a  force  could  be 
brought  together  to  repel  them? 
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The  writer  of  this  letter  was  Benjamin  Franklin.  On  10  September,  1944,  more  than  160 
years  later,  a  copy  of  this  quote  was  kept  on  the  desk  of  another  American  located  in 
England.  This  American  was  Lieutenant  General  Lewis  H.  Brereton  and  on  that  day  in 
September,  he  was  responsible  for  planning  Operation  MARKET  -  the  allied  airborne 
invasion  of  Holland  (11:122). 

Operation  MARKET-GARDEN  was  a  combined  airborne  and  land  based 
invasion.  The  Eirst  Allied  Airborne  Army  was  to  drop  in  Holland  and  hold  key  bridges 
along  the  route  to  and  across  the  Rhine.  This  was  Operation  MARKET.  The  British 
XXX  Corps  armored  unit  would  drive  up  a  narrow  corridor  of  advance,  relieving  the 
airborne  units  as  it  went,  until  it  crossed  the  Rhine.  This  was  Operation  GARDEN.  If 
successful,  it  could  have  brought  about  an  early  end  to  the  war.  However,  this  was  not  to 
be. 

MARKET-GARDEN  would  require  three  major  airdrops  of  troops  over  three 
days.  Once  on  the  ground,  the  airborne  units  would  require  airdrop  resupply.  While  it 
would  be  difficult  to  identify  any  one  element  that  led  to  the  failure  of  MARKET- 
GARDEN  as  being  decisive,  the  lack  of  precision  airdrop  capability  is  clearly  significant. 
History  records  abysmal  airdrop  accuracy.  British  airborne  troops  “watched  in  despair  as 
thirty-five  Stirling  bomber-cargo  planes  dropped  supplies  everywhere  but  on  the  [drop] 
zones.  Of  eighty-seven  tons  of  ammunition,  food  and  supplies  destined  for  the  men  of 
Arnhem,  only  twelve  tons  reached  the  troops.  The  remainder,  widely  scattered  to  the 
southwest,  fell  among  the  Germans”  (1 1 :376).  This  was  not  the  first  or  the  last  time  that 
airdrop  would  inadvertently  supply  the  enemy  rather  than  the  defenders. 
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Figure  1,  C-47’s  performing  low  altitude  airdrop  during  Operation  MARKET- 
GARDEN,  (Source:  http://www.qmfound.com/airborne2.gif) 

In  April  of  1972,  the  forces  of  North  Vietnam  launched  their  Easter  Offensive.  It 
was  an  effort  to  overrun  South  Vietnam  in  one  stroke.  The  key  to  South  Vietnam  was  the 
capitol  in  Saigon.  North  Vietnamese  forces  planned  to  launch  from  Cambodia  and  drive 
the  90  mile  distance  down  Highway  13  to  the  capitol.  On  this  highway,  approximately 
26  miles  from  the  Cambodian  border,  sat  the  city  of  An  Loc.  It  was  here  that  a  major 
battle  would  ensue  that  would  lead  to  a  two  month  long  siege.  By  1972,  the  majority  of 
American  ground  units  have  been  withdrawn.  The  Army  of  the  Republic  of  South 
Vietnam  (ARVN)  had  about  6,000  troops  in  An  Loc  to  defend  against  more  than  35,000 
North  Vietnamese  forces  (19).  It  would  fall  to  American  air  power  to  sustain  them. 
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Figure  2,  Route  of  attack  on  An  Loc  and  surrounding  area,  (Source: 
http://www,vnafmamn,com/Valiant_Anloc.html) 

Airdrop  crews  flying  in  support  of  the  forees  on  the  ground  at  An  Loo  faoed  a  lethal 

curtain  of  fire  including  .51  caliber,  37mm,  and  57mm  Anti  Aircraft  Artillery  (AAA)  (9). 

On  1 1  May,  the  first  SA-7  Strela,  Infrared  guided,  Man  Portable  Air  Defense 

(MANPAD)  weapon  was  fired  in  the  vicinity  of  An  Loc  (19).  Prior  to  this  time,  the  only 

technique  that  afforded  an  adequate  level  of  preoision  was  the  Low  Altitude  Parachute 

Extraction  System  (LAPES).  But  such  tactics  proved  to  be  suicidal  in  face  of  the  anti-air 

environment  around  An  Eoo.  The  374*^  Tactical  Airlift  Wing  had  an  operating 

detachment  at  Tan  Son  Nhut  Air  Base  and  was  tasked  with  the  airlift  mission  for  An  Loo 

(9).  They  developed  revolutionary  techniques  for  high  altitude  airdrop  called  Ground 

Radar  Aerial  Delivery  System  (GRADS)  and  Adverse  Weather  Aerial  Delivery  System 

(AW ADS).  These  techniques  allowed  for  improved  accuracy  for  airdrop  from  above 

12,000  feet  (2).  The  374'*’  also  developed  new  parachute  methodology.  They  devised  a 

method  for  airdrop  using  a  smaller  26  foot  diameter  “ring-slot”  high  velocity  parachute 
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than  the  standard  64  foot  diameter  G-12  paraehute  eanopy.  The  ring-slot  ehute  served 
not  to  deeelerate  the  load,  but  to  stabilize  it  at  it  fell.  Careful  paekaging  allowed  most 
types  of  loads  to  survive  the  landing  (9). 


Figure  3,  C-130  performing  LAPES  cargo  drop  during  the  siege  of  Khe  Sahn  in 
1968,  (Source:  http://www.qmfound,com/khe_sanhl.jpg) 


Figure  4,  26’  Ring  Slot  High  Velocity  Parachutes  in  flight,  (Source: 
http://www,pioneeraero,com/pop-ups/2-14-IMAGEl,htm) 


These  new  teehniques  allowed  the  defenders  of  An  Loo  to  hold  out  against  the 
overwhelming  odds  they  faoed.  The  Easter  Offensive  failed  and  South  Vietnam  survived 
for  another  three  years.  For  the  airdrop  foroes  involved,  the  final  tally  of  losses  were  15 
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aircrew  casualties,  numerous  wounded  airerew  members,  37  aireraft  damaged,  and  the 
loss  of  2  C-123  and  3  C-130E  aireraft  (19). 

The  techniques  developed  by  the  374*'^  in  the  support  of  An  Loo  would  last  far 
beyond  that  South  East  Asian  battlefield;  more  than  30  years  later,  versions  of  them  are 
still  in  use.  Airdrop  would  oontinue  to  play  a  key  in  military  operations  all  over  the 
world.  Operations  JUST  CAUSE,  PROVIDE  PROMISE,  ALLIED  EORCE, 

ENDURING  EREEDOM  and  others  would  see  airdrop  being  oalled  on  time  and  again. 

In  the  intervening  years,  the  threats  faoed  by  the  airorews  in  Vietnam  have  only 
intensified.  Reaching  an  adequate  level  of  precision  using  conventional  techniques  now 
plaees  aireraft  and  their  erews  at  unaeeeptably  high  risk. 

The  airdrop  mission  has  evolved  beyond  the  1970’s  era  solution.  Methods  are 
needed  to  operate  outside  the  Weapons  Engagement  Zone  (WEZ)  of  MANPADS  and 
AAA  while  reaching  totally  unprecedented  levels  of  preeision.  It  is  additionally  desirable 
that  airdrop  be  able  to  operate  at  an  offset  from  the  desired  PI.  Sueh  a  capability  will 
allow  covert  teams  to  be  resupplied  via  airdrop  without  their  position  being  highlighted 
by  overflight  of  the  drop  aireraft.  Eortunately,  threats  and  requirements  are  not  the  only 
thing  to  have  evolved  sinee  the  1970’s. 

Sinee  the  advent  of  the  Global  Positioning  System  (GPS)  in  the  1990’s  more  and 
more  military  systems  have  eome  to  rely  upon  the  navigation  teehnology.  The  Joint 
Direet  Attaek  Munition  (JDAM)  revolutionized  preeision  engagement  and  has  virtually 
beeome  a  household  word.  It  would  not  take  long  for  the  teehnology  that  made  the 
JDAM  possible  to  begin  transiting  to  the  airdrop  world.  The  stage  is  finally  set  for  the 
Joint  Preeision  Airdrop  System  -  JPADS. 


12 


Joint  Precision  Air  Drop  System  Overview 

In  traditional  airdrop,  the  aircrew  must  fly  the  aircraft  to  a  specific  point  in  the 
sky,  known  as  the  Computed  Air  Release  Point  (CARP).  The  CARP  is  calculated  using 
variables  such  as  payload  weight,  drop  altitude,  aircraft  velocity  vector,  wind  velocity 
vector,  and  location  of  the  intended  Point  of  Impact  (PI).  One  CARP  corresponds  to  one 
PI.  Miss  the  CARP  and  you  miss  the  PI.  Of  course,  hitting  the  CARP  does  not  guarantee 
that  you  will  hit  the  PI,  but  it  is  the  point  of  maximum  likelihood  given  the  quality  of  the 
data  input  into  the  calculations.  This  is  where  weather  sensitivities  become  important  to 
understand. 

JPADS  is  intended  to  revolutionize  how  airdrop  works.  This  is  about  more  than 
bringing  GPS  precision  to  CARP  calculations  though.  JPADS  is  a  Family  of  Systems 
that  allows  for  precision  airdrop  to  one  or  more  PI  from  medium  to  high  altitude  with  the 
option  of  significant  standoff  range  (i.e.,  without  the  need  to  fly  directly  over  the  PI  as  in 
traditional  airdrop).  These  capabilities  allow  for  significant  operational  flexibility.  For 
example  a  single  aircraft  could,  in  a  single  airdrop  pass,  drop  loads  to  different  Pis. 
Alternately,  one  or  more  aircraft  could  drop  loads  from  a  broad  Launch  Acceptability 
Region  (LAR)  to  hit  a  single  PI  from  various  launch  points. 
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Figure  5,  JPADS  Guided  Parachute  drop  capabilities  (6) 


These  capabilities  are  important  for  tactical  advantage  as  well  as  safety  since  JPADS 
allows  aircrews  to  drop  from  altitudes  and  standoff  ranges  which  are  safe  from  enemy 
surface-to-air  threats  and  terrain.  And  finally,  the  ability  to  drop  on  a  PI  without  direct 
over- flight  serves  to  further  protect  the  aircraft  and  crew  as  well  as  to  prevent 
highlighting  the  location  of  the  PI  and  the  airdrop  recipient.  Figure  6  shows  the  JPADS 
Systems  Architecture  Operational  View  (OV-1)  Diagram.  The  OV-1  is  a  graphic 
depicting  the  high-level  operational  concept  of  the  JPADS  architecture. 
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JPADS  Physical  Components. 


JPADS  consists  of  a  roll-on/roll-off  system  suite  for  the  aireraft,  a  mission 
planning  element,  and  a  variety  of  speeialized  Guided  Paraehute  systems.  The  Air  Force 
is  the  program  lead  for  developing  the  aircraft  systems  whieh  consists  of  the  JPADS  Fly- 
Away  Kit,  JPADS  Mission  Planner  (JPADS-MP)  software,  and  Global  Positioning 
System  (GPS)  Dropsondes.  The  US  Army  has  responsibility  for  the  development  of  the 
Guided  Paraehute  systems. 

JPADS  Fly-Away  Kit. 

The  JPADS  Fly-Away  Kit  is  self-contained  unit  designed  to  give  roll-on/roll-off 
JPADS  system  capability  for  an  airlift  platform  such  as  the  C-130  or  C-17.  The  kit 
contains  a  Precision  Air  Drop  System  (PADS)  software  configured  Panasonic  CF-29 
Toughbook  (also  known  as  the  PADS  Laptop  Computer  or  PLC),  a  Global  Positioning 
System  Retransmission  System  (GPS  -  RTS),  the  Advanced  PADS  Interfaee  Proeessor 
(APIP),  and  all  neeessary  eonneetions  for  the  system  and  the  aireraft.  Figure  7  shows  the 
JPADS  Fly-Away  Kit  in  its  stowed  and  unstowed  eonfiguration.  With  its  case,  the  Fly- 
Away  Kit  weighs  75  lbs.  The  Kit  is  developed  by  Planning  Systems  Inc,  Draper  Labs, 
and  the  Foreeast  Systems  Lab  of  the  NOAA  (6). 
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Figure  7.  JPADS  Fly-Away  Kit  components.  (1) 


JPADS  Guided  Parachute  Family. 

Presently,  several  system  types  are  under  consideration.  Among  these  are  the 
Affordable  Guided  Airdrop  System  (AGAS),  the  Screamer,  and  the  Sherpa.  Each  system 
differs  in  approach  to  the  guided  airdrop  problem  solution  as  well  as  in  overall 
performance  capabilities.  This  section  will  provide  background  on  each  system  and  how 
it  fits  into  the  JPADS  architecture. 

Affordable  Guided  Airdrop  System  (AGAS). 

Developed  in  joint  venture  by  Vertigo  Inc  and  Capewell,  AGAS  is  a  family  of 
systems  for  precision  airdrop  of  loads  from  200  to  10,000  pounds.  It  is  intended  to 
provide  high  accuracy  and  precision  at  low  cost  by  utilizing  off  the  shelf  parachutes  and 
rigging  components  and  is  essentially  a  strap-on  guidance  kit  for  the  standard  Container 
Delivery  System  (CDS).  The  AGAS  system  is  compatible  with  existing  inventory 
parachutes  such  as  the  G-12  and  the  26’  Ring  Slot  High  Velocity  Parachute.  The  heart  of 
the  system  is  the  Autonomous  Guidance  Unit  (AGU).  The  JPADS-MP  generates  a  Wind 
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Profile  which  is  used  to  calculate  a  wind  corrected  flight  trajectory.  This  trajectory  is 
passed  on  to  the  AGAS  AGU.  The  AGU  monitors  the  actual  flight  path  as  compared  to 
the  nominal  flight  path.  The  flight  path  is  then  adjusted  by  “slipping”  the  parachute 
control  risers.  The  Figure  below  displays  the  AGAS  mission  profile. 


Figure  8.  AGAS  flight  profile  with  Wind  Profile  corrections.  (Source:  http://www.vertigo- 

inc.com/agas) 

In  testing,  AGAS  has  proven  to  be  highly  accurate,  and  typically  has  the  smallest  Circular 
Error  Probable  (CEP)  of  the  candidate  systems.  However,  it  also  has  the  least  horizontal 
standoff  capability  among  the  candidates  at  approximately  5  km.  AGAS  loads  typically 
have  a  14-15  minute  Total  Time  Aloft.  The  following  figure  shows  example  AGAS  miss 
distances  and  their  associated  CEPs. 
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Figure  9.  Sample  drop  score  card  and  CEP  achieved  by  the  AGAS  Guided  Parachute  system. 
(Source:  http://www.vertigo-inc.com/agas/cep.jpg) 


Sherpa. 

The  Sherpa  guided  parachute  system  is  the  product  of  Mist  Mobility  Integrated 
System  Technology  (MMIST).  Sherpa  is  a  family  of  four  systems  with  load  capacities 
ranging  from  265  lbs  to  2200  lbs.  The  Sherpa  system  uses  a  large  Ram  Air  Parachute 
(RAP)  which  gives  the  system  a  significant  glide  range  and  maneuverability.  The  RAP 
affords  Sherpa  a  horizontal  standoff  range  of  up  to  20  km  from  a  drop  altitude  of  25,000 
feet.  A  unique  feature  of  Sherpa  is  the  option  to  provide  terminal  guidance  via  a  hand 
control  unit.  Otherwise,  the  Sherpa  uses  an  AGU  to  correct  for  windage  errors  with 
respect  to  the  preplanned  trajectory  calculated  by  the  JPADS-MP. 
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Figure  10,  MMIST  Sherpa  prepares  to  land,  (Source: 
http://www,mmist,ca/Sherpa,asp) 

The  following  chart  shows  a  series  of  Sherpa  2200  miss  distances  with  associated  100, 
200,  and  300  m  CEPs.  Sherpa  is  a  Commercial  Off-The-Shelf  (COTS)  system  already  in 
use  with  US  Marine  Corps. 


A  Sherpa  2200 
landing  points 

- 100m 

200m 

- 300m 

- 80%  Frequency 

of  last  20  drops 
(217m) 


Figure  11,  MMIST  Sherpa  sample  drop  score  card  and  CEP,  (Stoker,  2006) 
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Screamer. 


Strong  Enterprises  with  Robotek  Engineering  have  developed  the  Sereamer 
Precision  Cargo  Delivery  System.  The  Screamer  is  unique  in  that  it  uses  an  undersized 
Ram  Air  Drogue  (RAD)  rather  than  a  full  size  canopy.  The  RAD  serves  to  stabilize  and 
decelerate  the  payload  as  well  provide  steering  capability.  The  use  of  the  RAD  also 
allows  for  a  rapid  decent  from  altitude  and  improved  resistance  to  wind  effects. 

However,  due  to  its  small  size,  the  Screamer  RAD  is  incapable  of  slowing  the  payload 
down  for  landing.  This  is  accomplished  by  the  deployment  of  one  or  more  standard 
round,  unguided  parachutes  (typically  one  or  more  G-1 1  parachutes)  once  the  payload 
nears  the  surface.  Once  the  recovery  chute  is  deployed,  the  Screamer  is  considered  to  be 
ballistic,  which  is  to  say,  at  the  mercy  of  the  low-altitude  winds.  The  figures  below  show 
the  phases  of  Screamer  flight,  first  under  the  RAD  and  then  moments  after  the 
deployment  of  the  Recovery  Parachutes. 


First  USAF  JPADS 
Combat  Airdrop 


Screamer  2K:  31  Aug  06 


Figure  12,  Phases  of  Screamer  Flight,  (17, 10) 

As  with  the  other  systems.  Screamer  navigation  is  accomplished  via  an  AGE!  which 

compares  its  real-time  position  with  a  preplanned  trajectory  using  onboard  GPS.  The 
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Screamer  has  a  glide  ratio  of  2.6;  1 .  From  a  drop  altitude  of  25,000  feet,  it  has  a  standoff 
distance  of  about  7.2  miles  and  has  a  Total  Time  Aloft  of  approximately  8.6  minutes. 
The  following  Figure  show  example  miss  distances  for  the  Screamer  system  and  its 
associated  CEPs. 


Figure  13,  Scream  drop  score  card  and  CEP,  (10) 


Joint  Precision  Airdrop  System  -  Mission  Planner  (JPADS-MP). 

The  JPADS-MP  is  a  eombination  of  the  PLC  and  the  PADS  software.  Its  ultimate 
use  is  to  ensure  that  the  cargo  arrives  at  the  desired  PI.  Figure  14  shows  the  functional 
design  of  JPADS  as  well  as  the  data  flow  strueture. 
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Figure  14,  JPADS-MP  system  layout,  (1) 


The  box  at  the  center  of  the  Figure  shows  the  functions  performed  by  the  JPADS-MP. 
As  can  be  seen,  data  flows  into  the  JPADS-MP  from  the  GPS  dropsonde,  Combat  Track 
II  messages,  NIPRNET/SIPRNET,  and  from  user  input.  Data  flows  out  of  the  JPADS- 
MP  to  the  user  and  to  the  Guided  Airdrop  Systems.  A  more  comprehensive  picture 
comes  from  the  JPADS  SV-4  System  Functionality  Description  diagram  in  Figure  15. 


22 


Figure  15,  JPADS  SV-4,  Systems  Functionality  Description,  (6) 

These  two  diagrams  (Figures  14  and  15)  are  important  to  providing  the  answer  to 
the  first  Investigative  Question:  How  does  weather  affect  JPADS?  (i.e.,  how  does  the 
JPADS -MP  ingest  and  use  weather  data?  What  are  the  outputs?)  The  next  diagram  is  a 
simplified  version  of  the  SV-4.  All  non-weather  related  items  and  data  flows  have  been 
removed.  This  Systems  Engineering  product  answers  where  weather  enters  the  system 
and  how  that  data  flows  within  the  system  accomplish  its  functions. 
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Figure  16,  Weather  Data  Flows  within  JPADS, 

In  this  color  coded,  streamlined  version  of  the  SV-4,  the  data  flow  is  easy  to  follow. 
Looking  at  the  funetions  at  the  bottom  of  the  flow,  we  find  (as  one  would  expect)  the 
AGU  and  the  aireraft  navigation  system.  What  stands  out  as  the  eritieal  funetion  lies  one 
step  above  these;  Generate  CARP  -  the  Computed  Air  Release  Point,  a  funetion  of  the 
PLC.  The  flows  shown  here  are  the  primary  ones  available  to  operational  users.  There 
are  options  not  shown  here.  These  inelude  Pilot  Reports  (PIREPS),  Ballistie  Winds,  and 
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Climatology.  However,  it  turns  out  that  the  critieal  path  is  the  same  in  all  eases.  The 
PADS  Laptop  Computer  (PLC)  requires  at  least  two  inputs  to  ealeulate  the  CARP; 
mission  data  and  at  least  one  type  of  weather  data.  The  next  seetions  will  discus  these 
inputs  as  well  as  the  CARP  in  more  detail. 

Air  Force  Weather  Agency  (AFWA)  Weather  Forecasts. 

The  JPADS-MP  uses  4  Dimensional  forecast  models  in  order  to  generate  the  best 
CARP  solutions.  These  forecasts  are  called  4  Dimensional  as  they  include  x,  y,  and  z 
spatial  coordinates  as  well  as  a  temporal  coordinate.  AFWA  generates  forecasts  for  YPG 
in  three  levels  of  resolution:  5km,  15km,  and  45km.  In  4  Dimensional  forecast  models, 
resolution  refers  to  how  closely  spaced  the  data  weather  data  points  are  on  the  x/y  grid 
plane.  Thus,  5km  spacing  is  high  resolution,  containing  much  more  data  than  a 
resolution  with  45km  spacing  of  grid  points.  Obviously,  higher  resolution  means  a  larger 
data  file,  and  thus  greater  band-width  for  transmission  and  longer  download  times. 

The  5  and  15km  models  are  run  every  12  hours  and  the  45km  model  is  run  every 
6  hours.  One  model  run  simulates  24  hours  of  weather.  To  get  a  forecast  for  a  given 
time,  say  1200Z,  several  options  are  available.  The  5  and  15km  forecasts  are  initiated  at 
0600Z  and  1800Z.  To  get  a  forecast  valid  for  1200Z,  you  could  use  the  model  initiated 
on  the  current  day  at  0600  and  take  the  predicted  weather  conditions  6  hours  after  model 
start.  An  0600Z  start  time  plus  6  hours  lead  time  equals  a  valid  time  of  1200Z. 
Alternately,  the  1800Z  start  time  from  the  previous  day  with  an  18  hour  lead  time  also 
results  in  the  desired  1200Z  valid  time  on  the  current  day.  If  the  45km  resolution  is 
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considered,  five  more  forecasts  (three  from  the  current  day  and  two  from  the  previous) 
are  available  to  predict  the  conditions  at  1200Z. 

A  unique  aspect  of  weather  forecasting  is  that  models  are  initiated  “dry.”  This 
means  that  data  sueh  as  humidity,  dew  point,  and  pressure  are  fed  into  them,  but  not 
information  on  preeipitation  or  cloud  effects.  It  is  left  to  the  model’s  weather  physics  to 
generate  this  information.  This  results  in  a  certain  amount  of  spin-up  time  being  required 
by  the  model  before  it  begins  to  provide  realistic  forecast  results.  It  is  this  feature  whieh 
calls  into  question  how  model  lead  time  affects  forecast  accuracy. 

A  collection  of  these  forecasts  as  well  as  corresponding  weather  balloon 
soundings  covering  the  JPADS  ACTD  test  activity  from  20  June  2005  to  5  Deeember 
2006  has  been  provided  by  AFWA  Detachment  3  for  the  purposes  of  this  researeh.  A  list 
of  the  weather  balloons  specifically  used  for  this  analysis  is  included  in  Appendix  H. 

This  data  is  used  to  determine  whieh  forecasts  can  be  used  for  CARP  generation. 

Computed  Air  Release  Point  (CARP)  and  Launch  Acceptability  Region  (LAR). 

Certainly,  one  of  the  keys  to  preeise  airdrop  is  positioning  the  drop  aircraft  in  the 
proper  position  in  space  with  respeet  to  the  PI,  taking  into  aecount  the  variables  of 
aireraft  veloeity  as  well  as  the  wind  velocity  at  each  altitude  from  the  drop  level  down  to 
the  ground.  This  point  in  space  is  traditionally  known  as  the  Computed  Air  Release  Point 
(CARP).  One  of  the  ehief  funetions  of  the  JPADS-MP  is  ereating  a  highly  aecurate 
CARP.  This  is  accomplished  by  taking  into  account  aircraft  type,  altitude,  heading, 
airspeed,  position,  and  ramp-angle,  as  well  as  paraehute  type,  load  weight,  et  cetera.  To 
these  variables,  a  final  key  ingredient  is  added;  the  wind  profile.  While  the  payloads  are 
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guided,  they  are  unpowered  and  eannot  regain  kinetie  energy  onee  spent.  This  makes  a 
good  knowledge  of  the  air  mass  they  are  to  fly  through  eritieal  to  hitting  the  PI.  This 
CARP  is  then  input  (by  hand)  into  the  drop  aireraft  navigation  system.  Although  a 
preeise  formulation  of  the  CARP  is  not  as  eritieal  for  Guided  Paraehute  systems,  JPADS- 
MP  is  also  used  to  improve  the  aeeuraey  of  eheaper,  unguided  paraehute  systems  sueh  as 
the  High  Veloeity  Container  Delivery  System  (HV-CDS).  Thus,  operators  need  the  best 
possible  weather  estimate  to  ensure  aeeurate  airdrop. 

Sinee  the  Guided  Paraehute  systems  have  the  energy  to  fly  to  the  PI  from  a  large 
area,  the  JPADS-MP  caleulates  a  Launeh  Aeeeptability  Region  (LAR)  in  addition  to  the 
CARP.  The  LAR  is  an  elliptieal  region  whieh  represents  the  approximate  area  in  spaee 
from  whieh  a  Guided  Paraehute  system  eould  sueeessfully  reaeh  the  designated  PI  given 
the  weather  inputs  to  the  JPADS-MP.  Mathematieally,  it  is  the  solution  spaee  eontaining 
all  feasible  CARPS  for  the  Guided  Paraehute  systems  for  a  given  set  of  PI  eoordinates 
and  weather  inputs.  It  is  important  to  note  that  this  region  is  an  approximation  intended 
to  give  airerews  a  good  idea  of  the  system  limitations.  The  edge  of  the  LAR  should  not 
be  eonsidered  a  preeision  measurement  for  drop  purposes.  To  deal  with  this  and  other 
uneontrollable  variables  (sueh  as  the  true  weather  vs.  foreeast)  a  safety  faetor  of  1 1%  is 
subtraeted  from  the  LAR  ellipse.  The  safety  faetor  is  a  user  definable  option  within  the 
JPADS-MP. 

Figure  17  shows  a  eomparison  of  airdrop  missions  using  traditional,  non- JP ADS 
planning,  as  well  as  JPADS  planned  guided  and  unguided  drops.  Note  that  for  unguided 
drops,  only  one  CARP  is  available  to  hit  the  PI;  while  for  guided  drops  an  elliptieal  area 
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defined  by  an  Early,  Nominal,  and  Late  CARP  (the  Launeh  Aeeeptability  Region)  is 
sufficient  to  hit  the  PI. 


Tiadib'anal  Unguided  Airdrop  JPADS  Rsnned/Unguided  Airdrop  JPADS  Renned/Guided  Aintop 
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Figure  17,  Comparison  of  Flight  Profiles  for  traditional  vs,  guided  airdrop  options. 


Screamer  Recovery  Chute  Opening  Point  (OP). 

In  addition  to  the  CARP/LAR,  a  third  calculation  is  made  by  the  JPADS-MP  in 
support  of  the  Screamer  Guided  Parachute  system.  Unlike  the  AGAS  and  Sherpa,  the 
Screamer  requires  an  additional  AGU  command  beyond  the  CARP/LAR  computation. 
This  is  for  the  recovery  chute  Opening  Point  (OP),  also  called  the  pickle-point.  The 
correct  calculation  of  this  point  is  critical  to  hitting  the  PI  since  the  Screamer  payload  is 
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no  longer  guided  from  there.  Figure  18  demonstrates  the  partieular  ease  of  the  Sereamer 
system  and  its  sensitivity  to  correct  weather  forecasting.  The  conic  section  indicates  the 
volume  of  space  in  which  the  Screamer  system  has  sufficient  energy  to  maneuver  to 
reach  the  Opening  Point  (OP)  calculated  by  the  JPADS-MP.  Once  the  OP  is  reached,  the 
Screamer  payload  is  carried  to  the  ground  by  the  low-level  winds.  As  Figure  18  shows,  if 
these  winds  are  correctly  forecast.  Screamer  can  hit  the  PI  with  high  accuracy. 

Conversely,  poor  forecasting  leads  to  missing  the  PI. 
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Figure  18,  Screamer  Flight  Profile  and  weather  sensitivities. 
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GPS  Dropsonde. 

This  is  a  GPS  instrumented  unit  that  falls  under  paraehute  at  a  known  veloeity 
(typieally  70  fps).  The  dropsonde  is  released  from  an  aireraft  to  gather  a  sounding  of  the 
true  weather  in  elose  geographieal  and  temporal  proximity  to  the  planned  airdrop  loeation 
and  time.  The  weather  data  gathered  from  the  dropsonde  ean  be  integrated  in  flight  with 
the  preflight  mission  planning  foreeast  to  improve  the  preflight  planned  CARP/LAR  and 
OP. 

A  natural  assumption  is  that  updating  the  CARP/LAR  and  OP  generated  from  pre¬ 
mission  foreeasts  with  sampled  atmospherie  data  should  improve  their  estimates.  There 
are,  however,  potential  faults  in  this  assumption.  First,  dropsondes  take  time  -  both  to 
fall  and  for  their  data  to  be  assimilated  into  the  model.  In  order  to  ensure  adequate  time, 
eombat  taeties  eall  for  dropsondes  to  be  employed  no  later  than  ten  minutes  prior  to  the 
planned  airdrop.  Sinee  airdrop  missions  are  typieally  flown  at  approximately  150  mph, 
this  equates  to  a  minimum  differenee  of  more  than  25  miles  between  where  the 
atmosphere  was  sampled  and  where  the  aetual  drop  will  oeeur.  It  is  easily  eoneeivable 
that  a  dropsonde  eould  be  sampling  weather  on  one  side  of  a  mountain  ridge  and  the 
airdrop  take  plaee  on  the  other  side  in  eompletely  different  weather  eonditions. 
Additionally,  in  tests,  dropsonde  data  reeeption  beeomes  unreliable  at  low  altitude  - 
preeisely  the  time  when  aeeurate  information  is  most  oritieal. 
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Finally,  there  is  the  afore  mentioned  11%  factor  of  safety.  This  margin  of  safety 
is  essentially  energy  in  the  bank  for  the  on  board  guidance  system  to  use  should  it 
encounter  unexpected  weather  during  descent.  The  question  then  becomes,  is  a 
dropsonde  likely  to  ever  dictate  moving  the  LAR  more  than  11%,  particularly 
considering  the  other  limitations  of  dropsonde  employment?  Since  dropsondes  are  a 
consumable,  they  add  to  mission  cost  as  well  as  complexity. 

Initial  planning  for  this  study  called  for  an  analysis  of  these  questions.  However, 
this  was  de-scoped  from  the  thesis  after  consulting  with  Draper  Labs  concerning  the  LAR 
calculation  performed  by  the  JPADS-MP  version  used  in  this  research.  As  the  JPADS- 
MP  continues  to  evolve,  the  calculation  of  LAR  will  change  significantly  rendering  any 
work  done  valueless. 
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III.  Methodology 


Research  Strategy 

The  first  challenge  with  this  research  was  in  determining  what  was  meant  by 
weather  sensitivity  of  the  JPADS  system  and  how  to  measure  it.  The  intent  was  to  apply 
statistical  analysis  to  data,  but  to  what  data  and  how?  The  obvious  answer  to  the  first  part 
of  that  question  was  the  data  recorded  by  AFWA  Det  3  in  support  of  the  JPADS  ACTD. 
For  each  test,  a  record  was  kept  of  the  weather  balloons  launched  in  support  of  that  day’s 
missions  as  well  as  the  associated  valid  weather  forecasts.  This  data  covers  a  period  from 
20  June  2005  to  5  December  2006.  This  was  a  lot  of  data,  over  50  GB  worth.  The  next 
question  is  how  to  analyze  it.  AFWA  Det  3  has  already  begun  looking  at  a  direct, 
altitude  by  altitude,  comparison  between  weather  balloons  and  weather  forecasts.  This 
research  compared  the  forecast  wind  velocity  (heading  and  speed)  against  observed  wind 
velocity  for  the  three  different  forecast  resolutions  and  various  lead  times.  It  was  this 
initial  research  that  prompted  this  thesis. 

The  JPADS-MP  served  as  a  means  for  providing  an  apples-to-apples  analysis  of 
the  weather  input  options  within  the  context  of  the  Mission  Planner,  something  which 
would  be  very  useful  to  the  user  community.  This  is  possible  since  the  JPADS-MP  can 
perform  its  navigational  computations  from  either  the  forecasts  alone  or  the  weather 
balloons  alone.  This  capability  allows  for  the  creation  of  Northing  vs.  Easting  error 
comparisons,  much  like  the  Miss  Distance  charts  for  the  various  Guided  Parachute 
systems  shown  previously  in  Chapter  2.  To  execute  this,  a  standard  mission  scenario  is 
used  for  evaluating  all  input  data.  The  scenario  is  detailed  as  follows: 
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Mission  Name;  N  (for  due  North  Run-In  heading) 

Drop  Aireraft:  C-130 
Run-In;  360°  Magnetie 
Weather  Referenee  Point; 

YPG  Site  16  (Weather  Balloon  Launeh  Loeation) 

Lat;  N  33  19.800  Lon;  W  1 14  19.800 
Elevation;  1421  ftMSL 

Drop  Altitude;  17500  ft  AGE  Airspeed;  135  KIAS 
Magnetie  Variation;  12.346  W  (deg) 

Total  Ramp  Load(s);  1 
Eoads  To  Drop  This  Pass;  1 
Exit  Eoeation;  RAMP 
Stiek  Type;  Single 

Aireraft  Altimeter  Setting;  29.92  inehes  Hg 

Chute/System  Type;  Sereamer 

Total  Rigged  (All-up)  Weight  (lbs);  8000 

Plight  Station  (load  e.g.);  677 

Stiek  Position;  Eeft 

Glide  Safety  Paetor;  0.89 

PI;  YPG  JPADS  Center  PI 

PI  Coordinates;  N  33  19.612/W  114  22.226 
PI  Elevation;  1249  ft  MSP 
Ballistie  Chute  Type;  2  G1 1 

Steerable  Chute  Type;  850  Sq-Pt  (Sereamer  10k  System) 

This  N  mission  was  used  to  generate  CARP  and  OP  navigation  data  from  the  historieal 
weather  data.  These  were  then  grouped  by  resolution  for  analysis.  Analysis  was 
performed  in  Exeel,  Matlab,  and  JMP  6.  There  are  two  stages  to  the  analysis;  the  first 
eompares  the  three  resolutions,  and  the  seeond  eompares  lead  time.  The  eomparison 
variables  are  the  population  mean  and  varianee.  To  ensure  that  the  N  mission  was  not 
introdueing  error,  a  seeond  mission  was  tested  on  the  5  km  data  set.  This  S  mission 
differed  only  in  the  Run-In  heading  of  180°  Magnetie.  The  results  indieated  virtually  no 
differenee  from  the  CARPs  ealeulated  in  the  N  mission.  The  remainder  of  this  ehapter 
will  detail  how  the  N  Mission  was  entered  into  the  JPADS-MP  and  how  the  resulting  data 
was  eaptured  and  evaluated. 
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JPADS-MP  Operation  -  N  Mission 


The  JPADS-MP  is  developed  by  Draper  Labs  and  Planning  Systems  Ine  (PSI), 
and  a  eomplete  user’s  manual  is  available  from  them.  This  discussion  will  be  limited  to 
the  aspects  of  the  JPADS-MP  that  were  used  in  the  execution  of  this  research.  Appendix 
G  contains  a  sequence  of  figures  that  will  provide  the  reader  with  sufficient  familiarity 
with  the  JPADS-MP  Graphical  User  Interface  (GUI)  to  recreate  the  steps  taken  in  this 
research.  Upon  starting  the  JPADS-MP,  the  user  is  presented  with  the  main  GUI  page  as 
shown  in  following  figure: 


3  PADS  -  N  -  Release  Pass  1 


RIe  Launch  View  Help 


Figure  19.  JPADS-MP  main  GUI. 
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For  this  study,  the  coordinates  for  Site  16  at  YPG  are  used  for  the  weather 
forecast  reference  point  as  this  is  the  location  from  which  the  weather  balloons  were 
launched.  This  representative  mission  was  created  using  an  actual  test  point  from  the 
ACTD  program,  the  only  modification  being  a  change  in  the  Run-In  heading  to  a  cardinal 
direction.  As  a  result,  the  PI  is  set  as  being  the  JPADS  Center  PI  target  at  YPG,  as  used 
in  testing.  This  is  located  3.7  km  from  Site  16.  It  may  have  been  better  for  the  purposes 
of  this  analysis  to  have  set  Site  16  to  be  the  PI.  Unfortunately,  this  was  realized  too  late 
for  implementation.  Flowever,  any  error  incurred  by  this  is  believed  to  be  minimal  when 
considering  that  the  highest  weather  resolution  was  5  km. 


Figure  20,  Weather  GUI 

The  Weather  GUI,  shown  in  Figure  20,  is  where  most  of  the  work  in  this  research 
was  done.  The  next  step  is  to  acquire  weather  data.  This  is  done  by  selecting  one  of  the 
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options  under  the  Weather  Aequisition  section.  The  options  relevant  to  this  research  are; 
Dropsondes,  4D  Forecast,  Balloon,  and  Climatology. 

The  JPADS-MP  uses  the  4D  Forecasts  generated  by  AFWA.  These  come  in  a 
format  known  as  GRidded  Information  in  Binary  format  (GRIB)  files.  Once  these  are 
downloaded,  the  Browse  button  is  used  to  point  the  Weather  Source  GUI  to  the  location 
of  required  GRIB  files.  Once  the  appropriate  path  is  specified  in  the  “GRIB  Files 
Location”  field,  select  the  “Acquire  Forecast”  button.  This  will  read  the  weather  forecast 
into  the  JPADS-MP  Environmental  Data  folder. 


Figure  21,  Weather  GUI  with  4D  Forecast  loaded. 

In  Figure  21,  the  4D  Forecast  inventory  now  shows  an  increment  of  one  and  the 
Wind  File  Production  section  now  has  the  options  for  wind  file  generation  via  LAPS 
Forecast.  The  Local  Analysis  and  Prediction  System  (LAPS)  is  the  most  advanced 
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modeling  method  ineluded  within  the  JPADS-MP.  It  allows  eomplex  modeling  of  wind 
interaetion  with  terrain  features  sueh  as  how  wind  will  flow  over  or  around  terrain 
obstaeles.  Seleet  either  the  Best  Available  or  LAPS  Forecast-only  (available  under  Full 
Options)  to  begin  Wind  File  produetion. 
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Figure  22,  JPADS  Main  GUI  CARP  Solution  TAB  after  successful  CARP 

calculation. 


Selecting  Compute  CARP  will  now  automatically  open  the  CARP  solution  tab.  The 


CARP  section  shows  the  Latitude  and  Longitude  of  the  Early,  Nominal,  and  Late  CARPs 


which  also  define  the  boundaries  of  the  LAR.  In  order  to  collect  this  data,  an  Optical 
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Screen  Reader  tool  was  developed  by  Captain  Ryan  Eggert  of  the  Air  Foree  Researeh 
Laboratory  Advaneed  Arehiteeture  and  Integration  Braneh. 

The  Sereen  OCR  tool  reads  the  values  in  the  Early,  Nominal,  and  Late  CARP 
coordinate  boxes  and  eopies  them  to  a  text  file.  In  doing  so,  it  also  eonverts  them  from  a 
DDD  MM.mmm  format  to  a  DDD.dddddd  format.  The  eonversion  to  deeimal  degrees 
allows  for  easier  mathematieal  operations  later.  Additionally,  the  Screen  OCR  copies  the 
eoordinates  for  the  Sereamer  OP  from  its  memory  loeation  and  writes  it  to  the  same  text 
file. 


The  method  of  building  text  files  for  analysis  is  to  segregate  the  data  by  weather 
balloons.  The  Sereen  OCR  allows  for  a  new  file  to  be  opened  and  then  to  append 
subsequent  data  to  this  file.  First,  the  CARP/LAR/OP  is  ealeulated  for  a  weather  balloon. 
This  data  is  saved  to  a  new  file  bearing  the  date  and  time  of  the  balloon  launch  as  the  file 
name.  Next,  the  CARP/LAR/OP  is  ealeulated  for  eaeh  weather  foreeast  that  was  valid 
for  the  time  of  that  weather  balloon  launeh.  Eaeh  new  data  set  is  appended  to  the  text  file 
resulting  in  a  file  similar  to  the  one  shown  below: 


Figure  23,  Sample  text  file  record  of  CARP  and  OP  calculations  from  the  JPADS- 
MP  captured  by  the  Screen  OCR  program. 
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As  can  be  seen,  eaeh  line  represents  a  different  weather  input:  weather  balloon  on 
the  first  line,  followed  by  weather  forecasts  of  varying  resolution  and  initialization  time. 
The  eoordinates  of  the  CARPs  and  OP  are  to  the  left  of  the  metadata.  Capt  Eggert  also 
developed  a  CARP  Analysis  tool  to  generate  Northing  and  Easting  data  from  the  raw 
eoordinates  eaptured  by  the  Sereen  OCR. 

The  CARP  analysis  tool  functions  by  comparing  each  weather  forecast  to  the 
weather  balloon  data  in  line  one  of  the  text  file.  This  results  in  a  file  similar  to  the  one 
shown  in  Figure  24,  below: 


Figure  24,  Sample  text  file  containing  output  from  the  CARP  Analysis  Tool, 

In  this  file,  the  data  represents  error  in  the  foreeasting.  A  value  in  the  Nominal 
NS  column  of  -40.578024  means  that  particular  forecast  generated  a  Nominal  CARP 
coordinate  that  was  40.578024  m  South  of  the  eorrect  Nominal  CARP  eoordinate  as 
defined  by  the  Nominal  CARP  ealeulated  from  the  weather  balloon  (an  aetual  sampling 
of  the  atmosphere).  Also  note  that,  while  the  resolution  data  is  unchanged,  the 
initialization  time  has  been  replaeed  by  the  Lead-Time.  This  is  aecomplished  by  simply 
taking  the  difference  between  the  weather  balloon  launeh  time  and  the  foreeast 
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initialization  time.  The  data  from  each  weather  balloon  (and  its  corresponding  forecasts) 
is  saved  in  a  folder  named  for  the  day  the  balloons  were  launched  on. 

Once  all  the  data  has  been  run  through  the  JPADS-MP  and  the  final 
Northing/Easting  data  has  been  saved,  the  whole  lot  is  read  into  Microsoft  Excel.  Excel 
is  used  to  organize  the  data  into  continuous  columns  by  resolution  and  then  order  them 
according  to  Eead-Time.  The  first  order  of  business  was  to  determine  if  a  separate 
analysis  would  need  to  be  performed  on  the  Early,  Nominal,  and  Eate  CARPS.  However, 
comparing  scatter  diagrams  for  each  type  of  CARP  indicated  this  was  unnecessary  and 
that  the  Nominal  CARP  would  suffice  for  all. 

Each  resolution  is  then  entered  into  Matlab  to  test  for  Bivariate  Normality.  This 
test  is  taken  from  Walsh  and  Eynch’s  discussion  on  the  Multivariate  Normal  Distribution 
(16:2).  It  was  possible  to  code  the  test  they  describe  into  Matlab  to  produce  a  Goodness- 
of-Eit  test  for  scaled  distances  to  a  Chi-Squared  Distribution  with  n  degrees  of  freedom. 
These  are  then  fit  to  a  regression  model.  The  R  adjusted  for  the  fit  then  give  an 
indication  of  the  GoE,  where  linearity  correlates  to  normality.  The  Matlab  input  script 
and  function  are  included  in  Appendix  E  and  E,  respectively. 

Having  passed  this  test,  the  data  sets  are  then  entered  into  JMP  6  for  detailed 
analysis.  JMP  6  was  used  to  perform  Analysis  of  Variance  (ANOVA)  as  well  as 
Multivariate  Analysis.  This  was  first  performed  for  the  full  data  set  of  each  resolution  in 
order  to  characterize  each  and  determine  if  one  was  more  favorable  than  the  others  in 
terms  of  mean  (error)  and  variance.  Then,  each  set  was  subdivided  in  order  to  examine 
the  effect  of  Eead-Time  on  sample  mean  and  variance.  Eor  5  and  15  km  data,  there  was 
insufficient  data  for  anything  other  than  a  morning  vs.  afternoon  comparison.  The  45  km 
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data,  however,  was  suffieient  to  group  Lead-Times  into  seven  bins  of  three  hours  eaeh. 
The  Lead-Times  for  45  km  resolution  range  from  approximately  2  to  23  hours. 
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IV,  Results  and  Analysis 


First  Look 

Before  commencing  the  statistical  analysis,  the  first  objective  was  to  verify  that 
earlier  assumptions  made  in  setting  up  the  test  were  valid.  The  following  diagram  is  a 
scatter  plot  showing  CARP  data  generated  from  the  comparison  of  5km  weather  data  to 
their  corresponding  weather  balloon  derived  CARP.  The  points  on  the  graph  show  the 
error  in  the  forecast  based  CARP  with  respect  to  the  “true”  weather  balloon  based  CARP. 
The  first  check  was  to  ensure  that  it  would  not  be  necessary  to  test  the  Early,  Nominal, 
and  Late  CARPs  individually,  but  rather,  that  one  category  would  suffice  for  all.  This 
chart  shows  that  the  errors  for  each  type  of  CARP  are  perfectly  correlated  and  validates 
the  concept  of  analyzing  only  the  Nominal  CARP  as  a  representative  for  the  whole. 


Corrparison  of  Early,  Norrinal,  &  Late  CARP  Scatter  for  N  Mission,  SkmData 
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Figure  25,  Scatter  Plot  of  Northing  and  Easting  Errors  for  Early,  Nominal,  and 

Late  CARPs, 
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The  next  chart  compares  the  N  Mission  used  in  this  study  with  a  notional  S 
Mission.  The  only  difference  being  a  Run-In  heading  of  180°  magnetic  as  opposed  to 
360°  magnetic.  The  purpose  of  this  test  is  to  determine  if  the  aircraft  velocity  vector 
played  a  significant  role  in  the  observed  CARP  errors.  As  can  be  seen  below,  there  is 
excellent  correlation  between  the  N  and  S  Missions,  discounting  any  such  concerns. 


Corrparison  of  Northin^Easting  for  N&S  Missions,  5  km  Data 


«  N  Mssion  (Nominal) 
□  S  Mssion  (Nominal)] 


Figure  26,  Scatter  Plot  of  Northing  and  Easting  errors  of  Nominal  CARP 
comparing  results  from  N  and  S  Missions, 

Figure  27  shows  the  results  of  the  full  data  set.  The  upper  chart  shows  the  CARP  errors 

in  Northing  and  Easting  between  weather  balloon  and  weather  forecast  inputs;  the  lower 

chart  displays  the  same  errors  for  the  Screamer  OPs. 
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Northing/Easting  of  S/1&45l<mV\6c  Forecasts  vs.  V\6c  Balloon 
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Figure  27,  Scatter  Plot  of  Northing  and  Easting  errors  of  Nominal  CARPs  at  5, 15, 

and  45km  Resolution 


Northing/Easting  of  5/15/45km  Forecast  OPs  vs.  V\fc  Balloon  OPs 
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Figure  28,  Scatter  Plot  of  Northing  and  Easting  errors  of  OPs  at  5, 15,  and  45km 

Resolutions 
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Goodness  of  Fit  (GoF)  Testing  for  Bivariate  Normality 

The  next  assumption  to  eheck  is  that  of  Bivariate  Normal  distribution  of  the  data. 
As  mentioned  in  Chapter  3,  this  is  accomplished  by  fitting  a  line  to  a  comparison  of 
scaled  distances  to  a  Chi-Square  distribution.  As  can  be  seen  in  the  figure  below,  The 
CARP  error  data  is  a  good  fit  to  Bivariate  Normal.  However,  the  OP  data  is  strongly 
influenced  by  outliers  which,  when  included  in  the  line-fit  calculation,  cause  the  OP  data 
to  fail  the  GoF  test.  Exclusion  of  these  outliers  allows  for  fits  (shown  in  green  on  the 
charts)  with  R2  adjusted  in  the  range  of  0.98  -  0.99;  clearly  an  excellent  fit. 
Unfortunately,  using  historical  data,  there  is  no  way  to  account  for  the  cause  of  these 
outliers.  Therefore,  for  the  purposes  of  this  study,  they  will  not  be  removed. 


Figure  29,  Bivariate  Normal  Goodness-of-Fit  test  applied  to  CARP  data  (top  row) 
and  OP  data  (bottom  row)  at  5, 15,  and  45km  Resolution  (columns  1,2,  and  3 

respectively). 
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Since  this  establishes  Bivariate  Normal  as  a  good  distribution  to  deseribe  the  data, 
we  now  move  on  to  analyzing  the  data  in  that  light.  The  next  series  of  figures  will 
display  statistieal  data  neeessary  to  answer  whether  there  is  an  ideal  weather  foreeast 
resolution  for  ealeulating  the  CARP. 

Figure  30  shows  the  CARP  error  seatter  for  the  5,  15,  and  45km  resolution  data. 
The  green  lines  indieate  the  mean  value  for  Northing  and  Easting.  The  solid  red  line  is 
the  Least  Squares  regression  fit  and  the  broken  red  lines  indieate  the  95%  eonfidence 
interval  around  the  fit.  The  fit  is  indieative  of  the  eorrelation  between  Northing  and 
Easting.  The  aqua  line  and  shaded  region  is  the  95%  density  ellipse  for  the  data  set.  It  is 
worth  noting  that  in  both  the  5  and  15km  resolutions,  Flo  (there  is  no  eorrelation  between 
Northing  and  Easting)  eannot  be  rejeeted  at  a  =  0.05.  However,  for  the  45km  resolution. 
Ho  is  rejeeted  at  a  =  0.05.  This  ean  been  seen  in  Eigure  30,  as  the  95%  eonfidenee 
interval  for  the  45km  data  does  not  inelude  a  line  of  zero  slope. 


Figure  30,  Scatter  Plot  of  Northing  and  Easting  errors  with  mean  errors  (green 
lines),  correlation  (solid  red  line),  95%  confidence  interval  on  correlation  (dotted 
red  line),  and  95%  density  ellipse  displayed  for  5, 15,  and  45km  Resolutions. 
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This  positive  correlation  was  unexpected.  As  with  the  question  of  outlier  data  in 
the  OP  analysis,  there  is  no  clear  cut  answer  to  the  source  of  this  correlation.  It  does  call 
into  question  the  N  Mission  setup  as  a  possible  explanation.  Weather  data  (both  balloon 
and  forecast)  records  wind  direction  using  true  headings.  However,  aircrews  typically 
plan  using  magnetic  headings.  Since  the  N  Mission  borrowed  its  details  from  an  actual 
mission,  aircraft  Run-In  headings  were  entered  using  magnetic  headings.  The  magnetic 
variance  at  YPG  is  approximately  13°.  It  is  unknown  if  this  plays  a  role  in  the  observed 
correlation  or  not  as  there  was  insufficient  time  for  testing  after  the  discovery  of  the 
anomaly. 

What  is  clear  from  these  figures  is  the  effect  of  resolution  on  both  the  means  and 
variance  of  CARP  errors.  General  improvement  in  mean  error  is  seen  as  resolution 
decreases  from  5km  to  45km.  However,  the  finer  resolutions  (i.e.,  5  and  15km)  have 
lower  variance  than  does  the  45km  resolution.  Additionally,  all  three  resolutions  exhibit 
a  marked  Northing  error.  The  following  tables  provide  summary  statistics  for  Figure  30. 
The  sample  means  and  standard  deviations  are  contained  in  the  Correlation  Table. 
Complete  output  from  IMP  6  is  included  for  all  data  in  Appendix  B,  C,  and  D. 


47 


Linear  Fit 

5km  N/S  =  104.20727  +  0.0403265  5km  E/W 
1 5km  N/S  =  111 .32849  +  0.0493597  1 5km  E/W 
45km  N/S  =  72.951882  +  0.1476738  45km  E/W 

Summary  of  Fit  _ 
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Anaiysis  of  Variance 


Source 

DF 
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1 
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DF 

Sum  of  Squares 
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Model 

1 
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Error 
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Source 

DF 
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Model 

1 
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Error 
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0.0002 

Correiation 
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Mean 
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Number 
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Correlation 

SIgnIf.  Prob 

Number 

45km  E/W 

7.961625 

398.495 

0.188678 

0.0002 

375 

45km  N/S 

74.12761 

311.893 

Figure  31,  JMP  6  Statistical  Output  for  5, 15,  and  45km  Resolution  CARP  data. 

The  answer  to  the  question  of  resolution  appears  to  be  that  that  45km  data 


provides  the  lowest  mean  error,  but  the  greatest  varianee.  The  next  question  is  that  of  the 


effect  of  Lead-Time  on  CARP  error.  Figures  32  and  33  show  the  frequency  of  Lead- 


Times  for  each  data  set.  Recall  that  Lead-Time  is  the  delta  between  forecast  initialization 
and  the  planned  drop  time.  Since  the  historical  data  used  was  never  intended  for  this  type 
of  study,  it  presents  certain  difficulties  which  will  now  be  addressed. 
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Testing  at  YPG  typically  occurs  in  two  temporal  groups:  before  noon  and  after 
noon  local  time.  Since  the  5  and  15km  forecasts  are  only  generated  twice  per  day  (at 
0600Z  and  1800Z),  there  are  significantly  less  data  points  available  for  the  Lead-Time 
study  at  these  resolutions  than  for  the  45km  data  (which  is  generated  every  6  hours). 
Ideally,  there  would  be  at  least  thirty  data  points  for  each  hour  of  Lead-Team  to  allow  for 
a  complete  comparison;  unfortunately  that  is  not  the  case.  In  order  to  ensure  enough  data 
for  statistical  significance,  Lead-Times  must  be  grouped  together  in  “bins.”  Due  to  the 
general  paucity  of  data  at  the  5  and  15km  resolution,  these  were  lumped  into  two  bins  at 
the  natural  break  point  in  the  histogram.  This  compares  Short  Lead-Times  (8  to  17 
hours)  to  Long  Lead-Times  (17-24  hours). 


5km  Lead-Time  Histogram 
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Figure  32,  Lead-Time  Histogram  for  5  and  15km  Resolution  data. 

The  45km  data  is  more  extensive,  but  still  requires  grouping  for  best  results.  In 
this  case  there  are  seven  bins  containing  three  hours  of  data  each  with  a  range  from  2 
hours  of  Lead-Time  out  to  23  hours  -  a  much  more  complete  set  of  observations. 
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45km  Lead  Time  h-istogram 


Bin 


Figure  33,  Lead-Time  Histogram  for  45km  Resolution  Data. 

The  analysis  begins  as  before;  this  time  with  Northing  and  Easting  error  plots 
differentiated  by  Lead-Time  as  well  as  Resolution.  We  then  move  on  to  a  One-Way 
Layout  to  further  investigate  the  behavior  of  the  means  and  variances  as  Lead-Time  is 
adjusted. 
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Variable 

Mean 

Std  Dev 
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Signif.  Prob 

Number 

5km  E/W  LTB  1 
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83 
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Figure  34,  Scatter  Plot  of  Northing  and  Easting  errors  for  5  and  15km  Resolution 
sorted  by  Lead-Time  bins  as  well  as  the  associated  statistical  data. 
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Means  and  Std  Deviations 
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80 
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31.968 
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81 
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81 
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1 
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1 
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Positive  values  show  pairs  of  means  that  are  significantly  different. 


Figure  35.  Results  of  JMP  6  Oneway  Layout  Analysis  for  5km  Resolution  data. 
Northing  data  is  on  the  left  side  of  the  figure  and  Easting  is  on  the  right  for  ease  of 

comparison. 
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Figure  35  is  a  combination  of  the  JMP  6  output  for  the  5km  resolution  forecast 
data.  It  allows  for  a  side-by-side  comparison  of  the  means  and  variances  of  both  the 
Northing  and  Easting  data.  The  first  two  graphs  show  the  distance  errors  with  their 
associated  bin  number.  The  inner  set  of  blue  dashes  indicates  the  mean  of  that  set  and  its 
confidence  interval.  The  outer  set  of  blue  dashes  indicates  1  Standard  Deviation.  The 
red  rings  to  the  right  of  the  chart  are  a  visualization  tool  for  comparing  means.  When  this 
data  is  displayed  in  JMP  6,  selecting  one  ring  will  cause  it  to  be  highlighted  with  a  thick 
red  ring  (as  opposed  to  a  standard  thin,  black  ring).  Subsequently,  all  groups  whose 
means  are  not  significantly  different  change  from  black  rings  to  red  rings.  Groups  with 
significantly  different  means  become  gray.  This  test  indicates  that  there  is  no  significant 
difference  in  the  means  of  Lead-Time  bins  1  and  2  in  either  the  Northing  or  the  Easting 
data.  The  lower  portion  of  the  table  is  a  test  to  verify  that  the  variance  between  the  bins 
is  not  significant.  JMP  6  applies  five  different  methods  to  this  evaluation.  In  each  case, 
the  high  p-Value  indicates  failure  to  reject  Hq:  the  variances  between  bins  are  equal. 
Eigure  36  presents  the  same  analysis  for  the  15km  resolution  data  with  comparable 
results. 
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Means  and  Std  Deviations 


Level 

Number 

Northing  Mean 

Std  Dev 

Std  Err  Mean 

Lower  95% 

Upper  95% 

1 

80 

78.083 

276.278 

30.889 

16.600 

139.57 

2 

83 

145.552 

280.181 

30.754 

84.373 

206.73 

Level 

Number 

Easting  Mean 

Std  Dev 

Std  Err  Mean 

Lower  95% 

Upper  95% 

1 

80 

77.100 

374.591 

41.881 

-6.3 

160.46 

2 

83 

-30.151 

354.638 

38.927 

-107.6 

47.29 

Tests  that  the  Variances  are  Equal 


15km  LT  Bins 


400 


15km  LT  Bins 


Test 

F  Ratio 

0'Brien[.5] 

0.0168 

Brown-Forsythe 

0.1502 

Levene 

0.1432 

Bartlett 

0.0157 

F  Test  2-slded 

1.0285 

DFNum 

DFDen 

p-Value 

1 

161 

0.8970 

1 

161 

0.6989 

1 

161 

0.7056 

1 

0.9002 

82 

79 

0.9013 

Test 

F  Ratio 

0'Brien[.5] 

0.2584 

Brown-Forsythe 

0.1041 

Levene 

0.1018 

Bartlett 

0.2397 

F  Test  2-slded 

1.1157 

DFNum 

DFDen 

p-Value 

1 

161 

0.6119 

1 

161 

0.7473 

1 

161 

0.7501 

1 

0.6244 

79 

82 

0.6235 

Means  Comparisons  for  each  pair 

t  Alpha 
1.97481  0.05 


using  Student's  t 

t 

1.97481 


Alpha 

0.05 


Abs(Dlf)-LSD  2  1 

2  -85.305  -18.631 

1  -18.631  -86.889 


Abs{Dlf)-LSD  1  2 

1  -113.83  -5.55 

2  -5.55  -111.76 


Positive  values  show  pairs  of  means  that  are  significantly  different. 


Figure  36.  Results  of  JMP  6  Oneway  Layout  Analysis  for  15km  Resolution  data. 
Northing  data  is  on  the  left  side  of  the  figure  and  Easting  is  on  the  right  for  ease  of 

comparison. 
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Again,  means  and  variances  do  not  appear  to  vary  significantly  between  Lead- 
Times  at  15km  resolution.  We  move  on  now  to  the  45km  resolution  data.  Figures  37  and 
38  display  the  CARP  error  data  for  each  of  the  seven  Lead-Time  bins  of  the  45km  data. 

In  this  sequence  of  charts,  the  Easting  mean  remains  relatively  close  to  zero  with  the 
greatest  deviation  occurring  in  Lead-Time  bin  5,  which  represents  data  Lead-Time  of  15 
to  17  hours.  Of  more  interest  are  the  results  of  the  Northing  mean.  For  the  first  two  bins 
(2-8  hours  Lead-Time),  the  Northing  mean  is  very  close  to  zero.  The  precise  values  of 
the  means  are  highlighted  in  Figure  39.  The  trend  of  the  Northing  error  mean  is  generally 
worse  beyond  8  hours  of  Lead-Time. 


Figure  37,  Scatter  Plot  of  45km  Resolution  Northing  and  Easting  error  sorted  by 

bin,  Bins  1-4, 
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Figure  38,  Scatter  Plot  of  45km  Resolution  Northing  and  Easting  error  sorted  by 

bin,  Bins  5-7, 

Bin  5  is  again  the  location  of  largest  mean  error,  this  time  for  the  Northing  error. 
Additionally,  Bin  5  shows  the  greatest  level  of  potential  correlation  between  Northing 
and  Easting  error. 

The  next  series  of  tables  presents  the  Oneway  Layout  for  the  45km  resolution 
forecast  data.  There  are  some  differences  in  the  data  presented  here  due  to  the  addition 
levels  (i.e.,  bins)  available  for  comparison.  First  is  the  Tukey-Kramer  HSD  (Honestly 
Significant  Difference)  Test  in  addition  to  the  pair-wise  Student’s  t  Test. 
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Oneway  Analysis  of  45km  CARP  Error  By  45km  Lead-Time  Bins 


Means  and  Std  Deviations 


Level 

Number 

Northing  Mean 

Std  Dev 

Std  Err  Mean 

Lower  95% 

Upper  95% 

1 

48 

0.10 

286.139 

41.301 

-82.98 

83.19 

2 

51 

-6.76 

274.515 

38.440 

-83.97 

70.45 

3 

53 

44.09 

294.586 

40.464 

-37.10 

125.29 

4 

60 

49.32 

278.429 

35.945 

-22.60 

121.25 

5 

53 

162.64 

345.121 

47.406 

67.51 

257.77 

6 

57 

120.32 

308.757 

40.896 

38.39 

202.24 

7 

53 

138.93 

359.118 

49.329 

39.94 

237.91 

Level 

Number 

Easting  Mean 

Std  Dev 

Std  Err  Mean 

Lower  95% 

Upper  95% 

1 

48 

-15.632 

397.339 

57.351 

-131.0 

99.74 

2 

51 

-10.805 

359.155 

50.292 

-111.8 

90.21 

3 

53 

11.717 

378.416 

51.979 

-92.6 

116.02 

4 

60 

5.949 

363.862 

46.974 

-88.0 

99.94 

5 

53 

65.746 

413.535 

56.803 

-48.2 

179.73 

6 

57 

-12.650 

442.575 

58.620 

-130.1 

104.78 

7 

53 

10.295 

441.901 

60.700 

-111.5 

132.10 

Figure  39,  JMP  6  Oneway  Layout  Analysis  of  45km  Resolution  data,  Part  1, 
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Tests  that  the  Variances  are  Equal 


500-1 - 

400- 

>  ”  ■  " 

^  300--i - ; - ■ - ; - ^ - 

2  200  = 

100- 

o-' - 1 - 1 - 1 - 1 - 1 - 1 — 

1  2  3  4  5  6  7 

45km  LT  Bins 


500-1 - 

400 - ■ - i - = - = 

_  ■  ■ 

1  300- 

2  200  = 

100- 

o-' - 1 - 1 - 1 - 1 - 1 - 1 — 

1  2  3  4  5  6  7 

45km  LT  Bins 


Test 

F  Ratio 

DFNum 

DFDen 

Prob  >  F 

Test 

F  Ratio 

DFNum 

DFDen 

Prob  >  F 

0'Brien[.5] 

1 .2329 

6 

368 

0.2885 

0'Brien[.5] 

1 .0776 

6 

368 

0.3753 

Brown-Forsythe 

1.9188 

6 

368 

0.0768 

Brown-Forsythe 

0.8953 

6 

368 

0.4983 

Levene 

1.9868 

6 

368 

0.0667 

Levene 

1 .0443 

6 

368 

0.3961 

Bartiett 

1.1944 

6 

0.3057 

Bartiett 

0.7998 

6 

0.5699 

Weich  Anova  testing  Means  Equai, 

aiiowing  Std  Devs  Not  Equa 

F  Ratio 

DFNum 

DFDen 

Prob  >  F 

F  Ratio 

DFNum 

DFDen 

Prob  >  F 

2.4046 

6 

162.36 

0.2408 

6 

162.4 

0.9624  1 

Northing  Means  Comparisons  for  each  pair  using  Student's  t 


Levei 

-  Levei 

Difference 

Lower  CL 

Upper  CL 

5 

2 

169.3968 

50.542 

288.2513 

5 

1 

162.5347 

41.802 

283.2672 

7 

2 

145.6865 

26.832 

264.5410 

7 

1 

138.8244 

18.092 

259.5569 

6 

2 

127.0771 

10.286 

243.8686 

6 

1 

120.2150 

1.513 

238.9171 

5 

3 

118.5439 

0.838 

236.2500 

5 

4 

113.3135 

-0.908 

227.5349 

7 

3 

94.8336 

-22.872 

212.5398 

7 

4 

89.6032 

-24.618 

203.8247 

6 

3 

76.2242 

-39.398 

191.8468 

6 

4 

70.9938 

-41.079 

183.0670 

4 

2 

56.0833 

-59.321 

171.4878 

3 

2 

50.8529 

-68.002 

169.7074 

4 

1 

49.2212 

-68.117 

166.5589 

3 

1 

43.9908 

-76.742 

164.7233 

5 

6 

42.3197 

-73.303 

157.9424 

5 

7 

23.7103 

-93.996 

141.4164 

7 

6 

18.6095 

-97.013 

134.2321 

1 

2 

6.8621 

-114.990 

128.7144 

4 

3 

5.2304 

-108.991 

119.4518 

p-Vaiue 

Difference 

i  -  I  :  li 

i  .1  li 

i  1  i 

1  i 

0.0518 

0.1140 

0.1238 

0.1957 

0.2137 

0.3399 

0.4007 

0.4100 

0.4741 

0.4721 

0.6923 

0.7518 

0.9119 

0.9283 

1 _ _ ! _ 1 

1 _ !_^ _ ! _ 1 

1 

1 

1 

1  1  1  1 

1 _ ! _ H _ ! _ 1 

i  1  H 

1  1  1  1 

U _ ■ _ ! _ 1 

U _ ■ _ ! _ 1 

1 _ 1 _ ! _ 1 

n - i - 1 - i 

Figure  40,  JMP  6  Oneway  Layout  Analysis  of  45km  Resolution  data,  Part  2 


As  previously  mentioned,  the  rings  to  the  right  of  the  data  plots  are  a  visual  tool 
for  comparing  the  means  of  each  group  Based  on  the  Student’s  t  Test,  we  see  two 
groups  of  similar  means:  1,2,3,  and  4  in  the  first  group  and  5,6,  and  7  in  the  second.  A 
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green  line  has  been  added  to  the  plot  to  indieate  zero  mean  error.  The  blaek  line  on  the 
plot  displays  the  overall  mean  error.  As  ean  be  seen,  the  first  grouping  of  means  fall 
below  the  blaek  line  while  the  seeond  grouping  lies  above  it.  This  eorresponds  (as 
expeeted)  with  the  seatter  plot  data  indieating  generally  greater  error  as  Lead-Time 
inereases,  with  the  greatest  error  loeated  in  bin  5.  Note  the  data  table  for  the  Student’s  t 
eomparisons.  Cells  highlighted  in  red  indieate  a  p-Value  less  than  a  =  0.05  and  thus  a 
rejeetion  of  Hq:  no  signifieant  differenee  in  means  between  bins. 

The  Tukey-Kramer  HSD  Test  is  used  here  in  addition  to  the  Student’s  t  beeause 
the  number  of  observations  in  eaeh  bin  is  unequal.  In  sueh  conditions,  Tukey-Kramer  is 
intended  to  provide  a  conservative  test  for  difference  in  the  means.  In  this  case,  the 
Tukey-Kramer  result  differs  from  the  Student’s  t,  indicating  that  that  Ho  cannot  be 
rejected  under  that  test.  For  the  Easting  error  data,  both  Student’s  t  and  Tukey-Kramer 
fail  to  reject  Ho,  which  is  wholly  expected  from  a  visual  inspection  of  the  plot. 

In  both  cases,  bins  with  shorter  Lead-Times  display  smaller  Standard  Deviations, 
and  hence,  smaller  variance.  However,  the  testing  suite  for  equal  variance  fails  to  reject 
Ho:  no  significant  difference  in  variances  between  bins.  The  complete  IMP  6  report  is 
available  in  Appendices  B,  C,  and  D. 

Next  we  will  briefly  examine  the  Screamer  Opening  Point  (OP)  error  scatter  plots. 
The  improvement  in  overall  error  as  compared  to  the  CARP  calculations  is  immediately 
apparent.  These  diagrams  are  on  a  300x300m  plot  as  opposed  to  the  1250x1250m  plots 
for  CARP  errors.  The  radius  for  95%  density  ellipse  radius  is  about  150m  at  all 
resolutions.  Additionally,  it  is  not  possible  to  reject  Ho  (no  correlation  between  Northing 
and  Easting)  for  the  OPs  at  any  resolution.  The  mean  Northing  and  Easting  errors  remain 
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virtually  constant  across  all  three  resolutions,  as  can  be  seen  in  Figure  41.  It  is  worth 
recalling  that  this  data  includes  anomalous  outliers.  These  were  left  in  the  data  set  as 
there  is  no  method  available  for  determining  their  cause  (and  thus  justify  their  removal) 
from  the  historical  data. 


Variable 

Mean 

Std  Dev 

Correlation 

Signif.  Prob 

Number 

5km  OP  E/W 

15.65148 

60.22144 

0.165302 

0.0355 

162 

5km  OP  N/S 

-1.77951 

48.40718 

Variable 

Mean 

Std  Dev 

Correlation 

Signif.  Prob 

Number 

15km  OP  E/W 

16.04198 

58.63512 

0.113049 

0.1508 

163 

15km  OP  N/S 

0.552814 

50.58542 

Variable 

Mean 

Std  Dev 

Correlation 

Signif.  Prob 

Number 

45km  OP  E/W 

16.2265 

58.74914 

0.084928 

0.1006 

375 

45km  OP  N/S 

-0.13898 

58.04491 

Figure  41,  Scatter  Plot  of  Northing  and  Easting  errors  for  5, 15,  and  45km 
Resolution  OP  data  as  well  as  the  associated  statistical  data. 
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This  chapter  concludes  with  a  quick  look  at  Yuma  Proving  Ground  Site  16,  the 
location  where  this  weather  data  was  collected.  Figure  42  shows  the  locations  of  YPG 


Site  16  as  well  as  the  JPADS  Center  PI  drop  zone. 


Figure  42,  Satellite  view  of  Yuma  Proving  Ground  Site  16  and  JPADS  Center  PI, 

(Source:  Google  Earth) 

As  can  be  seen  in  the  figure,  the  average  wind  at  Site  16  follows  the  local  terrain. 


Wind  in  this  area  can  generally  be  expeeted  to  flow  up  the  valley  during  daylight  heating 


and  down  the  valley  during  night-time  cooling.  The  effect  is  shown  in  the  diagram 
below: 
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It  is  reasonable  to  believe  that  this  local  weather  effect  may  be  responsible  for  the  greater 
variability  observed  in  CARP  Easting  errors.  However,  certainty  would  require  testing  at 
other  locations.  Finally,  Figure  44  is  a  composite  image  showing  a  Bivarite  Normal 
distribution  generated  in  Matlab  from  the  45km  data  with  an  overlay  of  the  Eigen  vectors 
of  the  distribution.  The  Eigen  vectors  are  collinear  with  the  axes  of  the  ellipse.  The 
angle  between  the  Eigen  vectors  and  the  y-axis  is  approximately  18.7°.  While  this  value 
is  greater  than  the  magnetic  variance  for  the  area,  it  is  not  by  much.  Also,  it  is  aligned 
very  closely  to  the  direction  of  True  North,  as  can  been  in  Figure  44  below: 


Figure  44,  Overlay  of  Eigen  Vectors  on  45km  Resolution  derived  Bivariate  Normal 
distribution  along  with  the  approximate  relationship  to  True  and  Magnetic  North, 

This  argues  in  favor  of  the  idea  that  inputting  magnetic  heading  rather  than  true  may  be 

the  cause  of  the  apparent  correlation. 
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V.  Summary  and  Conclusions 


The  Bottom  Line 

This  study  set  out  to  determine  how  JPADS-MP  outputs  were  sensitive  to  weather 
inputs.  The  main  question  being:  are  there  ideal  weather  inputs  to  obtain  the  most 
accurate  outputs?  To  answer  this  question  required  gaining  a  detailed  understanding  of 
how  these  three  entities  (input-system-output)  interact.  Not  surprisingly,  some  questions 
remain.  However  there  are  clear,  useful  results  from  this  research: 

45km  resolution  weather  forecast  data  has  lower  mean  CARP  calculation  error 
than  does  5  or  15km  resolution  weather  forecast  data.  {More  accurate) 

45km  resolution  weather  forecast  data  has  higher  variance  in  CARP 
calculation  error  than  does  5  or  15km  resolution  forecast  data.  {Less  precise) 

2  to  8  hours  of  Lead-Time  offers  lowest  mean  CARP  calculation  error  with 
least  variance. 

Screamer  Opening  Point  CARP  calculation  errors  do  not  appear  to  depend 
upon  weather  forecast  resolution. 

For  the  JPADS-MP,  there  appears  to  be  no  advantage  to  using  higher 
resolution  weather  forecast  data.  High  resolution  weather  forecasts  have 
larger  file  sizes  requiring  greater  bandwidth  utilization  and  more  time  to 
acquire.  Further,  they  are  only  run  twice  daily  as  opposed  to  the  four  times 
daily  run  of  the  45km  resolution  data.  It  is  clear  that  within  the  context  of  the 
JPADS-MP  as-is,  the  lower  resolution  45km  weather  forecast  data  is  of 
greatest  value  to  the  user. 
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Future  Research 


This  study  actually  had  two  goals:  the  first,  to  attempt  to  answer  the  stated 
research  question;  the  seeond,  to  determine  if  further  researeh  on  this  particular  topic  was 
needed.  The  answer  to  the  seeond  question  is  elearly:  yes!  Now  that  there  is  a  better 
understanding  of  the  system,  the  best  way  to  achieve  definitive  results  is  to  eonduet  a 
Designed  Experiment  rather  than  attempting  to  make  do  with  historieal  data.  While  this 
is  a  time  intensive  process,  this  study  has  produeed  sufficiently  interesting  results  to 
warrant  the  effort.  Goals  for  a  follow-on  study  should  inelude  the  following: 

Gather  weather  balloon  and  weather  foreeast  data  from  several, 
geographieally  diverse  loeations. 

Gather  a  sufficient  number  of  weather  foreeasts  at  eaeh  resolution  to  allow  full 
eoverage  of  potential  Lead-Times  as  well  as  to  avoid  the  need  of  grouping 
Lead-Times  into  bins  (as  well  as  an  equal  number  of  observations  at  eaeh 
level). 

Determine  if  entering  navigation  data  using  magnetic  heading  versus  true 
heading  is  a  souree  of  error. 

Attempt  to  determine  eause  of  the  exeessive  Northing  errors. 

The  final  reeommendation  for  future  researeh  is  to  address  deficieneies  in  the 
method  for  aequiring  elimatology  data  for  JPADS.  JPADS-MP  uses  elimatology  data  as 
filler  whenever  data  neeessary  for  ealeulations  is  not  available  from  the  eurrent  weather 
input  data.  Additionally,  elimatology  may  be  used  as  a  last  resort  for  mission  planning 
should  no  other  souree  of  weather  input  be  available.  At  present,  elimatology  is  aequired 
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by  sending  a  request  to  the  Air  Foree  Combat  Climatology  Center  and  waiting  for  them 
to  provide  the  data.  This  ean  take  as  long  as  two  weeks.  AFIT’s  Engineering  Physies 
department  has  developed  an  excellent  software  alternative  known  as  the  Laser 
Environmental  Effects  Definition  and  Reference  (LEEDR),  which  could  allow  the 
JPADS-MP  user  to  have  the  ability  to  generate  accurate  climatology  data  on  the  same 
PADS  Laptop  Computer  that  they  use  for  the  JPADS-MP.  This  would  require  only 
minor  modifications  to  allow  the  two  programs  to  talk  to  each  other  (certainly  a  Systems 
Engineering  area  of  expertise!)  and  would  provide  immediate  and  tangible  improvement 
in  the  system  for  the  warfighter. 
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Appendix  A:  List  of  Acronyms 


AAA 

Anti-Aircraft  Artillery 

AGU 

Airborne  Guidance  Unit 

ACID 

Advaneed  Concept  Technology  Demonstration 

AFIT 

Air  Foree  Institute  of  Teenology 

AFRL 

Air  Foree  Researeh  Laboratory 

AFWA 

Air  Force  Weather  Ageney 

AGAS 

Affordable  Guided  Parachute  System 

AGL 

Above  Ground  Level 

ANOVA 

Analysis  of  Variance 

APIP 

Advanced  PADS  Interfaee  Proeessor 

AWADS 

Adverse  Weather  Aerial  Delivery  System 

CARP 

Computed  Air  Release  Point 

CDS 

Container  Delivery  System 

CEP 

Cireular  Error  Probable 

DoD 

Department  of  Defense 

GoF 

Goodness  of  Pit 

GPS 

Global  Positioning  System 

GRP-RTS 

GPS  Re-Transmit  Kit 

GRADS 

Ground  Radar  Air  Drop  System 

GRIB 

Gridded  Information  in  Binary  Pormat 

GUI 

Graphieal  User  Interface 

HVCDS 

High  Veloeity  Container  Delivery  System 
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JPADS 


JPADS-MP 

JDAM 

LAPPS 

LAPS 

LAR 

LEEDR 

MM5 

MMIST 

NOAA 

OP 

PADS 

PI 

PIREP 

PEC 

PSI 

RAD 

WEZ 

WRE 

YPG 


Joint  Precision  Air  Drop  System 
JPADS-Mission  Planner 
Joint  Direct  Attack  Munition 
Eow  Altitude  Parachute  Extraction  System 
Eocal  Analysis  and  Prediction  System 
Eaunch  Acceptability  Region 

Easer  Environmental  Effects  Definition  and  Reference 
Mesoscale  Model  5 

Mist  Mobility  Integrated  Systems  Technology 

National  Oceanic  and  Atmospheric  Administration 

Opening  Point 

Precision  Air  Drop  System 

Point  of  Impact 

Pilot  Report 

PADS  Eaptop  Computer 
Planning  System,  Inc. 

Ram  Air  Drogue 
Weapons  Engagement  Zone 
Weather  Research  and  Eorecasting 
Yuma  Proving  Ground 
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Appendix  B:  5km  JMP  6  Analysis  Output 
Bivariate  Fit  of  Exp  Chi-Sqr  By  Distance  (Nominal  CARP) 


- Linear  Fit 


Linear  Fit 

Exp  Chi-Sqr  =  0.2131093  +  0.7523327  Distance 


Summary  of  Fit 


RSquare 

0.97758 

RSquare  Adj 

0.97744 

Root  Mean  Square  Error 

0.29726 

Mean  of  Response 

1 .995723 

Observations  (or  Sum  Wgts) 

162 

Lack  Of  Fit 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Lack  Of  Fit 

158 

14.137851 

0.089480 

655.7718 

Pure  Error 

2 

0.000273 

0.000136 

Prob  >  F 

Totai  Error 

160 

14.138123 

0.0015 

Max  RSq 

1.0000 

Analysis  of  Variance 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

616.46791 

616.468 

6976.517 

Error 

160 

14.13812 

0.088 

Prob  >  F 

C.  Totai 

161 

630.60603 

<.0001 

Parameter  Estimates 


Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

0.2131093 

0.031638 

6.74 

<.0001 

Distance 

0.7523327 

0.009007 

83.53 

<.0001 
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Bivariate  Fit  of  Exp  Chi-Sqr  By  Distance  (Nominai  OP) 


- Linear  Fit 

- Linear  Fit 


Linear  Fit 

Exp  Chi-Sqr  =  0.8372231  +  0.4746891  Distance 

Summary  of  Fit 

RSquare  0.820416 

RSquareAdj  0.819294 

Root  Mean  Square  Error  0.841303 

Mean  of  Response  1 .995723 

Observations  (or  Sum  Wgts)  1 62 


Anaiysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

517.35942 

517.359 

730.9491 

Error 

160 

113.24661 

0.708 

Prob  >  F 

C.  Totai 

161 

630.60603 

<.0001 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio  Prob>|t| 

intercept 

0.8372231 

0.078773 

10.63  <.0001 

Distance 

0.4746891 

0.017558 

27.04  <.0001 

Linear  Fit 

Exp  Chi-Sqr  =  0.060006  +  0.9201943  Distance 

Summary  of  Fit 

RSquare 

0.995629 

RSquare  Adj 

0.995601 

Root  Mean  Square  Error 

0.10011 

Mean  of  Response 

1.746387 

Observations  (or  Sum  Wgts) 

156 

Anaiysis  of  Variance 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

351.59272 

351.593 

35082.19 

Error 

154 

1 .54338 

0.010 

Prob  >  F 

C.  Totai 

155 

353.13610 

<.0001 

Parameter  Estimates 


Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

0.060006 

0.012054 

4.98 

<.0001 

Distance 

0.9201943 

0.004913 

187.30 

<.0001 
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Bivariate  Fit  of  5km  N/S  By  5km  E/W  (Nominai  CARP) 


Linear  Fit 

Bivariate  Normal  Eiiipse  P=0.95(l 


Linear  Fit 

5km  N/S  =  104.20727  +  0.0403265  5km  E/W 


Summary  of  Fit 


RSquare  0.002742 

RSquare  Adj  -0.00349 

Root  Mean  Square  Error  286.6814 

Mean  of  Response  105.679 

Observations  (or  Sum  Wgts)  1 62 


Anaiysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

36162 

36162.2 

0.4400 

Error 

160 

13149793 

82186.2 

Prob  >  F 

C.  Totai 

161 

13185955 

0.5081 

Parameter  Estimates 


Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

104.20727 

22.63283 

4.60 

<.0001 

5km  E/W 

0.0403265 

0.060794 

0.66 

0.5081 

Correiation 

Variabie 

Mean 

Std  Dev 

Correiation 

Signif.  Prob 

Number 

5km  E/W 

5km  N/S 

36.49546 

105.679 

371.6411 

286.1824 

0.052369 

0.5081 

162 
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Linear  Fit 

Bivariate  Normai  Eilipse  P=0.95(l 


Linear  Fit 

5km  N/S  LTB  1  =  58.347947  +  0.1315542  5km  E/W  LTB  1 

Summary  of  Fit 

RSquare 
RSquare  Adj 
Root  Mean  Square  Error 
Mean  of  Response 
Observations  (or  Sum  Wgts) 


Analysis  of  Variance 


Source 

DF 

Sum  of  Squares  Mean  Square 

F  Ratio 

Modei 

1 

196129.5 

196129 

2.4428 

Error 

78 

6262642.5 

80290 

Prob  >  F 

C.  Totai 

79 

6458772.0 

0.1221 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

58.347947 

32.56149 

1.79 

0.0770 

5km  E/W  LTB  1 

0.1315542 

0.084171 

1.56 

0.1221 

Correlation 

Variabie 

Mean 

Std  Dev 

Correiation 

Signif.  Prob 

Number 

5km  E/W  LTB  1 

89.39787 

378.7505 

0.174259 

0.1221 

80 

5km  N/S  LTB  1 

70.10861 

285.9311 

0.030366 

0.017935 

283.3554 

70.10861 

80 
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Linear  Fit 

Bivariate  Normai  Eilipse  P=0.95(l 


Linear  Fit 

5km  N/S  LTB  2  =  139.95308  -  0.0283622  5km  E/W  LTB  2 

Summary  of  Fit 

RSquare 
RSquare  Adj 
Root  Mean  Square  Error 
Mean  of  Response 
Observations  (or  Sum  Wgts) 


Analysis  of  Variance 


Source 

DF 

Sum  of  Squares  Mean  Square 

F  Ratio 

Modei 

1 

8415.6 

8415.6 

0.1033 

Error 

80 

6518795.5 

81484.9 

Prob  >  F 

C.  Totai 

81 

6527211.2 

0.7488 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

139.95308 

31.55152 

4.44 

<.0001 

5km  E/W  LTB  2 

-0.028362 

0.088254 

-0.32 

0.7488 

Correlation 

Variabie 

Mean 

Std  Dev 

Correiation 

Signif.  Prob 

5km  E/W  LTB  2 

-15.1166 

359.3861 

-0.03591 

0.7488 

5km  N/S  LTB  2 

140.3818 

283.8712 

0.001289 

-0.01119 

285.4557 

140.3818 

82 


Number 

82 
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Means  and  Std  Deviations 


Level 

Number 

Mean 

Std  Dev 

Std  Err  Mean 

Lower  95% 

Upper  95% 

1 

80 

70.109 

285.931 

31.968 

6.478 

133.74 

2 

82 

140.382 

283.871 

31.348 

78.008 

202.76 

Means  Comparisons 

Comparisons  for  each  pair  using  Student's  t 

t  Alpha 

1.97490  0.05 

Abs{Dif)-LSD  2  1 

2  -87.868  -18.142 

1  -18.142  -88.960 

Positive  values  show  pairs  of  means  that  are  significantly  different. 

Tests  that  the  Variances  are  Equal 

400  q - 

>  300- - ^ ^ - 

°200- 
co  100- 

0=1 - ^ - 

1  2 

5km  LT  Bins 


Level 

Count 

Std  Dev 

MeanAbsDif  to  Mean 

MeanAbsDif  to  Median 

1 

80 

285.9311 

229.3837 

229.3326 

2 

82 

283.8712 

236.1517 

236.1347 

Test 

F  Ratio 

DFNum 

DFDen 

p-Value 

0'Brien[.5] 

0.0048 

1 

160 

0.9448 

Brown-Forsythe 

0.0712 

1 

160 

0.7899 

Levene 

0.0706 

1 

160 

0.7908 

Bartlett 

0.0042 

1 

0.9486 

F  Test  2-sided 

1.0146 

79 

81 

0.9479 

Welch  Anova  testing  Means  Equal,  allowing  Std  Devs  Not  Equal 

F  Ratio  DFNum  DFDen  Prob  >  F 
2.4634  1  159.84  0.1185 

t  Test 
1 .5695 
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Oneway  Analysis  of  5km  E/W  By  5km  LT  Bins  (Nominal  CARP) 


Means  and  Std  Deviations 


Level 

Number 

Mean 

Std  Dev 

Std  Err  Mean 

Lower  95% 

Upper  95% 

1 

80 

89.398 

378.751 

42.346 

5.11 

173.68 

2 

82 

-15.117 

359.386 

39.688 

-94.08 

63.85 

Means  Comparisons 

Comparisons  for  each  pair  using  Student's  t 

t  Alpha 

1.97490  0.05 

Abs{Dif)-LSD  1  2 

1  -115.25  -10.03 

2  -10.03  -113.83 

Positive  values  show  pairs  of  means  that  are  significantly  different. 


Tests  that  the  Variances  are  Equal 


Level 

Count 

Std  Dev 

MeanAbsDif  to  Mean 

MeanAbsDif  to  Median 

1 

80 

378.7505 

305.3941 

305.3941 

2 

82 

359.3861 

294.4267 

292.4927 

Test 

F  Ratio 

DFNum 

DFDen 

p-Value 

0'Brien[.5] 

0.2339 

1 

160 

0.6293 

Brown-Forsythe 

0.1441 

1 

160 

0.7047 

Levene 

0.1079 

1 

160 

0.7430 

Bartlett 

0.2189 

1 

0.6399 

F  Test  2-sided 

1.1107 

79 

81 

0.6392 

Welch  Anova  testing  Means  Equal,  allowing  Std  Devs  Not  Equal 

F  Ratio  DFNum  DFDen  Prob  >  F 
3.2430  1  159.05  0.0736 

t  Test 
1.8008 
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Bivariate  Fit  of  5km  OP  N/S  By  5km  OP  E/W  (Nominai  OP) 


Linear  Fit 

Bivariate  Normal  Ellipse  P=0.95(l 


Linear  Fit 

5km  OP  N/S  =  -3.859172  +  0.1328734  5km  OP  E/W 


Summary  of  Fit 


RSquare 

0.027325 

RSquare  Adj 

0.021246 

Root  Mean  Square  Error 

47.8902 

Mean  of  Response 

-1.77951 

Observations  (or  Sum  Wgts) 

162 

Anaiysis  of  Variance 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Model 

1 

10308.70 

10308.7 

4.4948 

Error 

160 

366955.36 

2293.5 

Prob  >  F 

C.  Total 

161 

377264.06 

0.0355 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

Intercept 

-3.859172 

3.888374 

-0.99 

0.3225 

5km  OP  E/W 

0.1328734 

0.062673 

2.12 

0.0355 

Correiation 

Variable 

Mean 

Std  Dev 

Correlation 

Signif.  Prob 

5km  OP  E/W 

15.65148 

60.22144 

0.165302 

0.0355 

5km  OP  N/S 

-1.77951 

48.40718 

Number 

162 
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Appendix  C:  15km  JMP  6  Analysis  Output 


- Linear  Fit 


Linear  Fit 

Exp  Chi-Sqr  =  0.2023848  +  0.7264998  Distance 


Summary  of  Fit 


RSquare 

0.98215 

RSquare  Adj 

0.98204 

Root  Mean  Square  Error 

0.265244 

Mean  of  Response 

1 .995749 

Observations  (or  Sum  Wgts) 

163 

Lack  Of  Fit 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Lack  Of  Fit 

160 

11.327014 

0.070794 

891.8347 

Pure  Error 

1 

0.000079 

0.000079 

Prob  >  F 

Totai  Error 

161 

11.327093 

0.0267 

Max  RSq 

1.0000 

Anaiysis  of  Variance 

Source 

DF  Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

623.26073 

623.261 

8858.846 

Error 

161 

11.32709 

0.070 

Prob  >  F 

C.  Totai 

162 

634.58782 

<.0001 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio  Prob>|t| 

intercept 

0.2023848 

0.02819 

7.18  <.0001 

Distance 

0.7264998 

0.007719 

94.12  <.0001 
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Bivariate  Fit  of  Exp  Chi-Sqr  By  Distance  (Nominal  OP) 


- Linear  Fit 

- Linear  Fit 


Linear  Fit 

Exp  Chi-Sqr  =  0.8487173  +  0.5074542  Distance 


Summary  of  Fit 


RSquare 

0.793745 

RSquare  Adj 

0.792464 

Root  Mean  Square  Error 

0.901644 

Mean  of  Response 

1 .995749 

Observations  (or  Sum  Wgts) 

163 

Lack  Of  Fit 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Lack  Of  Fit 

160 

130.88671 

0.818042 

5728.384 

Pure  Error 

1 

0.00014 

0.000143 

Prob  >  F 

Totai  Error 

161 

130.88685 

0.0105 

Max  RSq 

1.0000 

Analysis  of  Variance 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

503.70097 

503.701 

619.5875 

Error 

161 

130.88685 

0.813 

Prob  >  F 

C.  Totai 

162 

634.58782 

<.0001 

Parameter  Estimates 


Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

0.8487173 

0.084327 

10.06 

<.0001 

Distance 

0.5074542 

0.020387 

24.89 

<.0001 

Linear  Fit 

Exp  Chi-Sqr  =  -0.067938  +  1 .0699483  Distance 

Summary  of  Fit 


RSquare  0.993236 

RSquareAdj  0.993193 

Root  Mean  Square  Error  0.12467 

Mean  of  Response  1 .747532 

Observations  (or  Sum  Wgts)  1 57 
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Lack  Of  Fit 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Lack  Of  Fit 

154 

2.4089644 

0.015643 

109.5384 

Pure  Error 

1 

0.0001428 

0.000143 

Prob  >  F 

Total  Error 

155 

2.4091072 

0.0760 

Max  RSq 

1.0000 

Analysis  of  Variance 

Source 

DF  Sum  of  Squares 

Mean  Square 

F  Ratio 

Model 

1 

353.77827 

353.778 

22761.81 

Error 

155 

2.40911 

0.016 

Prob  >  F 

C.  Total 

156 

356.18737 

<.0001 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

Intercept 

-0.067938 

0.015614 

-4.35 

<.0001 

Distance 

1.0699483 

0.007092 

150.87 

<.0001 
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Bivariate  Fit  of  15km  N/S  By  15km  E/W  (Nominai  CARP) 


Linear  Fit 

Bivariate  Normai  Eiiipse  P=0.95(i 


Linear  Fit 

15km  N/S  =  111.32849  +  0.0493597  15km  E/W 


Summary  of  Fit 


RSquare 

0.00421 1 

RSquare  Adj 

-0.00197 

Root  Mean  Square  Error 

279.7441 

Mean  of  Response 

112.4385 

Observations  (or  Sum  Wgts) 

163 

Lack  Of  Fit 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Lack  Of  Fit 

160 

12599340 

78745.9 

Pure  Error 

1 

0 

0.0 

Prob  >  F 

Totai  Error 

161 

12599340 

Max  RSq 

1.0000 

Anaiysis  of  Variance 

Source 

DF  Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

53276 

53275.6 

0.6808 

Error 

161 

12599340 

78256.8 

Prob  >  F 

C.  Totai 

162 

12652615 

0.4105 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio  Prob>|t| 

intercept 

111.32849 

21.95251 

5.07  <.0001 

15km  E/W 

0.0493597 

0.059823 

0.83  0.4105 

Correiation 

Variabie 

Mean 

Std  Dev  Correiation  Signif.  Prob  Number 

15km  E/W 

22.48796 

367.396 

0.064889  0.4105  163 

15km  N/S 

112.4385 

279.4684 
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Bivariate  Fit  of  15km  N/S  LTB  1  By  15km  E/W  LTB  1  (Nominal  CARP) 


Linear  Fit 

Bivariate  Normai  Eiiipse  P=0.95(i 


Linear  Fit 

15km  N/S  LTB  1  =  67.491734  +  0.1373675  15km  E/W  LTB  1 


Summary  of  Fit 


RSquare 

RSquare  Adj 

Root  Mean  Square  Error 

Mean  of  Response 

Observations  (or  Sum  Wgts) 


0.034689 

0.022313 

273.1787 

78.08283 

80 


Analysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

209175.4 

209175 

2.8030 

Error 

78 

5820874.9 

74627 

Prob  >  F 

C.  Totai 

79 

6030050.3 

0.0981 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

67.491734 

31.19057 

2.16 

0.0335 

15km  E/W  LTB  1 

0.1373675 

0.082049 

1.67 

0.0981 

Correlation 

Variabie 
15km  E/W  LTB  1 
15km  N/S  LTB  1 


Mean 

Std  Dev 

Correiation 

Signif.  Prob 

Number 

77.10047 

374.5915 

0.186249 

0.0981 

80 

78.08283 

276.2784 

80 


Bivariate  Fit  of  15km  N/S  LTB  2  By  15km  E/W  LTB  2  (Nominal  CARP) 


Linear  Fit 

Bivariate  Normal  Ellipse  P=0.95(l 


Linear  Fit 

15km  N/S  LTB  2  =  145.1 1784  -  0.0144124  15km  E/W  LTB  2 

Summary  of  Fit 

RSquare 
RSquare  Adj 
Root  Mean  Square  Error 
Mean  of  Response 
Observations  (or  Sum  Wgts) 


Analysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Model 

1 

2142.2 

2142.2 

0.0270 

Error 

81 

6434986.0 

79444.3 

Prob  >  F 

C.  Total 

82 

6437128.2 

0.8700 

Parameter  Estimates 


Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

Intercept 

145.11784 

31.05097 

4.67 

<.0001 

15km  E/W  LTB  2 

-0.014412 

0.087769 

-0.16 

0.8700 

0.000333 

-0.01201 

281.8586 

145.5524 

83 


Correlation 

Variable 

15km  E/W  LTB  2 
15km  N/S  LTB  2 


Mean 

Std  Dev 

Correlation 

Signif.  Prob 

Number 

-30.1506 

354.6377 

-0.01824 

0.8700 

83 

145.5524 

280.1813 
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Oneway  Analysis  of  15km  N/S  By  15km  LT  Bins  (Nominal  CARP) 


Means  and  Std  Deviations 


Level 

Number 

Mean 

Std  Dev 

Std  Err  Mean 

Lower  95% 

Upper  95% 

1 

80 

78.083 

276.278 

30.889 

16.600 

139.57 

2 

83 

145.552 

280.181 

30.754 

84.373 

206.73 

Means  Comparisons 

Comparisons  for  each  pair  using  Student’s  t 

t  Alpha 

1.97481  0.05 

Abs{Dif)-LSD  2  1 

2  -85.305  -18.631 

1  -18.631  -86.889 

Positive  values  show  pairs  of  means  that  are  significantly  different. 


Tests  that  the  Variances  are  Equai 


Level 

Count 

Std  Dev 

MeanAbsDif  to  Mean 

MeanAbsDif  to  Median 

1 

80 

276.2784 

222.6416 

221.5912 

2 

83 

280.1813 

232.0265 

231.3428 

Test 

F  Ratio 

DFNum 

DFDen 

p-Value 

0’Brien[.5] 

0.0168 

1 

161 

0.8970 

Brown-Forsythe 

0.1502 

1 

161 

0.6989 

Levene 

0.1432 

1 

161 

0.7056 

Bartlett 

0.0157 

1 

0.9002 

F  Test  2-sided 

1 .0285 

82 

79 

0.9013 

Welch  Anova  testing  Means  Equal,  allowing  Std  Devs  Not  Equal 

F  Ratio  DFNum  DFDen  Prob  >  F 
2.3960  1  160.91  0.1236 

t  Test 
1.5479 
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Oneway  Analysis  of  15km  E/W  By  15km  LT  Bins  (Nominal  CARP) 


1000- 


500- 


LU 

E 


-500- 


-1000- 


u 

1  2 

15km  LT  Bins 


Student's  t 
0.05 


Means  and  Std  Deviations 


Level 

Number 

Mean 

Std  Dev 

Std  Err  Mean 

Lower  95% 

Upper  95% 

1 

80 

77.100 

374.591 

41.881 

-6.3 

160.46 

2 

83 

-30.151 

354.638 

38.927 

-107.6 

47.29 

Means  Comparisons 

Comparisons  for  each  pair  using  Student’s  t 

t  Alpha 

1.97481  0.05 

Abs{Dif)-LSD  1  2 

1  -113.83  -5.55 

2  -5.55  -111.76 

Positive  values  show  pairs  of  means  that  are  significantly  different. 

Tests  that  the  Variances  are  Equai 

4nn-[—  ■ 

> 

<D 

°200- 

■o 


15km  LT  Bins 


Level 

Count 

Std  Dev 

MeanAbsDif  to  Mean 

MeanAbsDif  to  Median 

1 

80 

374.5915 

299.8042 

299.8042 

2 

83 

354.6377 

289.1837 

289.0077 

Test 

F  Ratio 

DFNum 

DFDen 

p-Value 

0’Brien[.5] 

0.2584 

1 

161 

0.6119 

Brown-Forsythe 

0.1041 

1 

161 

0.7473 

Levene 

0.1018 

1 

161 

0.7501 

Bartlett 

0.2397 

1 

0.6244 

F  Test  2-sided 

1.1157 

79 

82 

0.6235 

Welch  Anova  testing  Means  Equal,  allowing  Std  Devs  Not  Equal 

F  Ratio  DFNum  DFDen  Prob  >  F 
3.5185  1  159.66  0.0625 

t  Test 
1 .8758 
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Linear  Fit 

Bivariate  Normal  Eiiipse  P=0.95(l 


Linear  Fit 

15km  OP  N/S  =  -1.011746  +  0.0975291  15km  OP  E/W 


Summary  of  Fit 


RSquare 

0.01278 

Rsquare  Adj 

0.006648 

Root  Mean  Square  Error 

50.41698 

Mean  of  Response 

0.552814 

Observations  (or  Sum  Wgts) 

163 

Lack  Of  Fit 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Lack  Of  Fit 

160 

409241.40 

2557.76 

Pure  Error 

1 

0.00 

0.00 

Prob  >  F 

Totai  Error 

161 

409241.40 

Max  RSq 

1.0000 

Anaiysis  of  Variance 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

5297.85 

5297.85 

2.0842 

Error 

161 

409241.40 

2541.87 

Prob  >  F 

C.  Totai 

162 

414539.25 

0.1508 

Parameter  Estimates 

Term 

Estimate  Std  Error 

t  Ratio 

Prob>|t| 

intercept 

-1.011746  4.094969 

-0.25 

0.8052 

15km  OP  E/W 

0.0975291  0.067556 

1.44 

0.1508 

Correiation 

Variabie 

Mean  Std  Dev 

Correiation 

Signif.  Prob 

15km  OP  E/W 

16.04198  58.63512 

0.113049 

0.1508 

15km  OP  N/S 

0.552814  50.58542 

Number 

163 
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Appendix  D:  45km  JMP  6  Analysis  Output 


Bivariate  Fit  of  Exp  Chi-Sqr  By  Distance  (Nominal  CARP) 


- Linear  Fit 


Linear  Fit 

Exp  Chi-Sqr  =  -0.043957  +  0.8737868  Distance 


Summary  of  Fit 


RSquare 

0.996655 

RSquare  Adj 

0.996646 

Root  Mean  Square  Error 

0.115216 

Mean  of  Response 

1.998153 

Observations  (or  Sum  Wgts) 

375 

Lack  Of  Fit 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Lack  Of  Fit 

371 

4.9511034 

0.013345 

64.8420 

Pure  Error 

2 

0.0004116 

0.000206 

Prob  >  F 

Totai  Error 

373 

4.9515150 

0.0153 

Max  RSq 

1.0000 

Analysis  of  Variance 

Source 

DF  Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

1475.3287 

1475.33 

111137.2 

Error 

373 

4.9515 

0.013275 

Prob  >  F 

C.  Totai 

374 

1480.2802 

0.0000 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio  Prob>|t| 

intercept 

-0.043957 

0.008539 

-5.15  <.0001 

Distance 

0.8737868 

0.002621 

333.37  0.0000 
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Bivariate  Fit  of  Exp  Chi-Sqr  By  Distance  (Nominal  OP) 


- Linear  Fit 

- Linear  Fit 


Linear  Fit 

Exp  Chi-Sqr  =  0.6325086  +  0.6302233  Distance 

Summary  of  Fit 


RSquare  0.851123 

RSquare  Adj  0.850724 

Root  Mean  Square  Error  0.768656 

Mean  of  Response  1 .9981 53 

Observations  (or  Sum  Wgts)  375 


F  Ratio 
16924.85 
Prob  >  F 
<.0001 
Max  RSq 
1.0000 


F  Ratio 
2132.419 
Prob  >  F 
<.0001 


Parameter  Estimates 


Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

0.6325086 

0.049499 

12.78 

<.0001 

Distance 

0.6302233 

0.013648 

46.18 

<.0001 

Linear  Fit 

Exp  Chi-Sqr  =  -0.056533  +  1.0449162  Distance 

Summary  of  Fit 

RSquare 
RSquare  Adj 
Root  Mean  Square  Error 
Mean  of  Response 
Observations  (or  Sum  Wgts) 


0.988978 

0.988947 

0.145563 

1.65801 

354 


Lack  Of  Fit 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Lack  Of  Fit 

369 

220.37995 

0.597236 

Pure  Error 

4 

0.00014 

0.000035 

Totai  Error 

373 

220.38009 

Analysis  of  Variance 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Modei 

1 

1259.9001 

1259.90 

Error 

373 

220.3801 

0.59 

C.  Totai 

374 

1480.2802 
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Lack  Of  Fit 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Lack  Of  Fit 

348 

7.4582215 

0.021432 

607.3446 

Pure  Error 

4 

0.0001412 

0.000035 

Prob  >  F 

Total  Error 

352 

7.4583626 

<.0001 

Max  RSq 

1.0000 

Analysis  of  Variance 

Source 

DF  Sum  of  Squares 

Mean  Square 

F  Ratio 

Model 

1 

669.21491 

669.215 

31583.83 

Error 

352 

7.45836 

0.021 

Prob  >  F 

C.  Total 

353 

676.67327 

0.0000 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

Intercept 

-0.056533 

0.012366 

-4.57 

<.0001 

Distance 

1.0449162 

0.00588 

177.72 

0.0000 
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Bivariate  Fit  of  45km  N/S  By  45km  E/W  (Nominai  CARP) 


Linear  Fit 

Bivariate  Normal  Ellipse  P=0.95(l 


Linear  Fit 

45km  N/S  =  72.951882  +  0.1476738  45km  E/W 


Summary  of  Fit 

RSquare 

RSquare  Adj 

Root  Mean  Square  Error 

Mean  of  Response 

Observations  (or  Sum  Wgts) 


0.035599 

0.033014 

306.7014 

74.12761 

375 


Lack  Of  Fit 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Lack  Of  Fit 

371 

35086535 

94572.9 

Pure  Error 

2 

0 

0.0 

Prob  >  F 

Total  Error 

373 

35086535 

Max  RSq 

1.0000 

Anaiysis  of  Variance 

Source 

DF  Sum  of  Squares 

Mean  Square 

F  Ratio 

Model 

1 

1295163 

1295163 

13.7687 

Error 

373 

35086535 

94066 

Prob  >  F 

C.  Total 

374 

36381698 

0.0002 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio  Prob>|t| 

Intercept 

72.951882 

15.84116 

4.61  <.0001 

45km  E/W 

0.1476738 

0.039798 

3.71  0.0002 

Correiation 

Variable 

Mean 

Std  Dev 

Correlation  Signif.  Prob  Number 

45km  E/W 

7.961625 

398.495 

0.188678  0.0002  375 

45km  N/S 

74.12761 

311.893 
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Bivariate  Fit  of  45km  N/S  LTB  1  By  45km  E/W  LTB  1  (Nominai  CARP) 


Linear  Fit 

Bivariate  Normal  Ellipse  P=0.95(l 


Linear  Fit 

45km  N/S  LTB  1  =  2.0917632  +  0.1272016  45km  E/W  LTB  1 


Summary  of  Fit 


RSquare 

RSquare  Adj 

Root  Mean  Square  Error 

Mean  of  Response 

Observations  (or  Sum  Wgts) 


0.0312 

0.010139 

284.6842 

0.103331 

48 


Anaiysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean 

Square 

F  Ratio 

Model 

1 

120061.9 

120062 

1.4814 

Error 

46 

3728074.7 

81045 

Prob  >  F 

C.  Total 

47 

3848136.7 

0.2298 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

Intercept 

2.0917632 

41.12309 

0.05 

0.9597 

45km  E/W  LTB  1 

0.1272016 

0.104509 

1.22 

0.2298 

Correiation 

Variable 

45km  E/W  LTB  1 
45km  N/S  LTB  1 


Mean 

Std  Dev 

Correlation 

Signif.  Prob 

Number 

-15.6321 

397.339 

0.176635 

0.2298 

48 

0.103331 

286.1385 
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Bivariate  Fit  of  45km  N/S  LTB  2  By  45km  E/W  LTB  2  (Nominai  CARP) 


Linear  Fit 

Bivariate  Normal  Ellipse  P=0.95(l 


Linear  Fit 

45km  N/S  LTB  2  =  -4.657001  +  0.1945243  45km  E/W  LTB  2 


Summary  of  Fit 


RSquare 

RSquare  Adj 

Root  Mean  Square  Error 

Mean  of  Response 

Observations  (or  Sum  Wgts) 


0.064771 

0.045685 

268.1706 

-6.75877 

51 


Anaiysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Model 

1 

244051.7 

244052 

3.3936 

Error 

49 

3523859.1 

71915 

Prob  >  F 

C.  Total 

50 

3767910.9 

0.0715 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

Intercept 

-4.657001 

37.56873 

-0.12 

0.9019 

45km  E/W  LTB  2 

0.1945243 

0.105595 

1.84 

0.0715 

Correiation 

Variable 

45km  E/W  LTB  2 
45km  N/S  LTB  2 


Mean 

Std  Dev 

Correlation 

Signif.  Prob 

Number 

-10.8047 

359.1552 

0.254502 

0.0715 

51 

-6.75877 

274.5145 
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Bivariate  Fit  of  45km  N/S  LTB  3  By  45km  E/W  LTB  3  (Nominai  CARP) 


Linear  Fit 

Bivariate  Normal  Eiiipse  P=0.95(l 


Linear  Fit 

45km  N/S  LTB  3  =  43.369608  +  0.0618356  45km  E/W  LTB  3 


Summary  of  Fit 


RSquare 

RSquare  Adj 

Root  Mean  Square  Error 

Mean  of  Response 

Observations  (or  Sum  Wgts) 


0.006309 

-0.01317 

296.5201 

44.09412 

53 


Anaiysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

28472.0 

28472.0 

0.3238 

Error 

51 

4484132.2 

87924.2 

Prob  >  F 

C.  Totai 

52 

4512604.3 

0.5718 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

43.369608 

40.75006 

1.06 

0.2922 

45km  E/W  LTB  3 

0.0618356 

0.108663 

0.57 

0.5718 

Correiation 

Variabie 

45km  E/W  LTB  3 
45km  N/S  LTB  3 


Mean 

Std  Dev 

Correiation 

Signif.  Prob 

Number 

11.71682 

378.4157 

0.079432 

0.5718 

53 

44.09412 

294.5859 
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Bivariate  Fit  of  45km  N/S  LTB  4  By  45km  E/W  LTB  4  (Nominai  CARP) 


Linear  Fit 

Bivariate  Normal  Eiiipse  P=0.95(l 


Linear  Fit 

45km  N/S  LTB  4  =  48.626014  +  0.1 174245  45km  E/W  LTB  4 


Summary  of  Fit 


RSquare 

RSquare  Adj 

Root  Mean  Square  Error 

Mean  of  Response 

Observations  (or  Sum  Wgts) 


0.023548 

0.006713 

277.4933 

49.32453 

60 


Anaiysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean 

Square 

F  Ratio 

Modei 

1 

107706.7 

107707 

1 .3987 

Error 

58 

4466146.2 

77003 

Prob  >  F 

C.  Totai 

59 

4573852.8 

0.2418 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

48.626014 

35.8291 

1.36 

0.1800 

45km  E/W  LTB  4 

0.1174245 

0.099286 

1.18 

0.2418 

Correiation 

Variabie 

45km  E/W  LTB  4 
45km  N/S  LTB  4 


Mean 

Std  Dev 

Correiation 

Signif.  Prob 

Number 

5.948603 

363.862 

0.153455 

0.2418 

60 

49.32453 

278.4294 
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Bivariate  Fit  of  45km  N/S  LTB  5  By  45km  E/W  LTB  5  (Nominai  CARP) 


Linear  Fit 

Bivariate  Normal  Eiiipse  P=0.95(l 


Linear  Fit 

45km  N/S  LTB  5  =  146.66501  +  0.242952  45km  E/W  LTB  5 

Summary  of  Fit 

RSquare 
RSquare  Adj 
Root  Mean  Square  Error 
Mean  of  Response 
Observations  (or  Sum  Wgts) 


Anaiysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

524891.8 

524892 

4.7223 

Error 

51 

5668765.7 

111152 

Prob  >  F 

C.  Totai 

52 

6193657.4 

0.0344 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

146.66501 

46.38147 

3.16 

0.0026 

45km  E/W  LTB  5 

0.242952 

0.111801 

2.17 

0.0344 

Correiation 


Variabie 

Mean 

Std  Dev 

Correiation 

Signif.  Prob 

Number 

45km  E/W  LTB  5 

65.74553 

413.5354 

0.291113 

0.0344 

53 

45km  N/S  LTB  5 

162.638 

345.1214 

0.084747 

0.066801 

333.3951 

162.638 

53 
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Linear  Fit 

Bivariate  Normal  Eiiipse  P=0.95(l 


Linear  Fit 

45km  N/S  LTB  6  =  121.5945  +  0.1008843  45km  E/W  LTB  6 


Summary  of  Fit 


RSquare 

RSquare  Adj 

Root  Mean  Square  Error 

Mean  of  Response 

Observations  (or  Sum  Wgts) 


0.020912 

0.00311 

308.2769 

120.3183 

57 


Analysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean 

Square 

F  Ratio 

Modei 

1 

111637.0 

111637 

1.1747 

Error 

55 

5226906.7 

95035 

Prob  >  F 

C.  Totai 

56 

5338543.7 

0.2832 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

121.5945 

40.84925 

2.98 

0.0043 

45km  E/W  LTB  6 

0.1008843 

0.093081 

1.08 

0.2832 

Correlation 

Variabie 

45km  E/W  LTB  6 
45km  N/S  LTB  6 


Mean 

Std  Dev 

Correiation 

Signif.  Prob 

Number 

-12.6501 

442.5745 

0.144608 

0.2832 

57 

120.3183 

308.7574 
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Linear  Fit 

Bivariate  Normal  Eiiipse  P=0.95(l 


Linear  Fit 

45km  N/S  LTB  7  =  137.36229  +  0.1520644  45km  E/W  LTB  7 


Summary  of  Fit 


RSquare 

RSquare  Adj 

Root  Mean  Square  Error 

Mean  of  Response 

Observations  (or  Sum  Wgts) 


0.035013 

0.016092 

356.2166 

138.9278 

53 


Analysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

234806.0 

234806 

1.8505 

Error 

51 

6471403.7 

126890 

Prob  >  F 

C.  Totai 

52 

6706209.7 

0.1797 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

intercept 

137.36229 

48.94365 

2.81 

0.0071 

45km  E/W  LTB  7 

0.1520644 

0.111786 

1.36 

0.1797 

Correlation 

Variabie 

45km  E/W  LTB  7 
45km  N/S  LTB  7 


Mean 

Std  Dev 

Correiation 

Signif.  Prob 

Number 

10.2948 

441.9014 

0.187118 

0.1797 

53 

138.9278 

359.1178 

95 


Means  and  Std  Deviations 


Level 

Number 

Mean 

Std  Dev 

Std  Err  Mean 

Lower  95% 

Upper  95% 

1 

48 

0.10 

286.139 

41.301 

-82.98 

83.19 

2 

51 

-6.76 

274.515 

38.440 

-83.97 

70.45 

3 

53 

44.09 

294.586 

40.464 

-37.10 

125.29 

4 

60 

49.32 

278.429 

35.945 

-22.60 

121.25 

5 

53 

162.64 

345.121 

47.406 

67.51 

257.77 

6 

57 

120.32 

308.757 

40.896 

38.39 

202.24 

7 

53 

138.93 

359.118 

49.329 

39.94 

237.91 

Means  Comparisons 

Comparisons  for  each  pair  using  i 

Student's  t 

t 

Alpha 

1 .96643 

0.05 

Abs{Dlf)-LSD 

5 

7  6 

4  3 

1  2 

5 

-117.71 

-94.00  -73.30 

-0.91  0.84 

41.80  50.54 

7 

-94.00 

-117.71  -97.01 

-24.62  -22.87 

18.09  26.83 

6 

-73.30 

-97.01  -113.50 

-41.08  -39.40 

1.51  10.29 

4 

-0.91 

-24.62  -41.08 

-110.63  -108.99 

-68.12  -59.32 

3 

0.84 

-22.87  -39.40 

-108.99  -117.71 

-76.74  -68.00 

1 

41.80 

18.09  1.51 

-68.12  -76.74 

-123.68  -114.99 

2 

50.54 

26.83  10.29 

-59.32  -68.00 

-114.99  -119.99 

Positive  values  show  pairs  of  means  that  are  significantly  different. 

Level 

Mean 

5 

A 

162.6380 

7 

A 

B 

138.9278 

6 

A 

B 

120.3183 

4 

A 

B 

C 

49.3245 

3 

B 

C 

44.0941 

1 

C 

0.1033 

2 

C 

-6.7588 

Levels  not  connected  by  same  letter  are  significantly  different. 


Level 

-  Level 

Difference 

Lower  CL 

Upper  CL 

p-Value 

5 

2 

169.3968 

50.542 

288.2513 

0.0053 

5 

1 

162.5347 

41.802 

283.2672 

0.0085 

7 

2 

145.6865 

26.832 

264.5410 

0.0164 

7 

1 

138.8244 

18.092 

259.5569 

0.0243 

6 

2 

127.0771 

10.286 

243.8686 

0.0330 

6 

1 

120.2150 

1.513 

238.9171 

0.0472 

5 

3 

118.5439 

0.838 

236.2500 

0.0484 

- r 

— r~ 

~r~ 

~i — 

1 — 
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Level 

-  Level 

Difference 

Lower  CL 

Upper  CL 

p-Value 

5 

4 

113.3135 

-0.908 

227.5349 

0.0518 

7 

3 

94.8336 

-22.872 

212.5398 

0.1140 

7 

4 

89.6032 

-24.618 

203.8247 

0.1238 

6 

3 

76.2242 

-39.398 

191.8468 

0.1957 

6 

4 

70.9938 

-41.079 

183.0670 

0.2137 

4 

2 

56.0833 

-59.321 

171.4878 

0.3399 

3 

2 

50.8529 

-68.002 

169.7074 

0.4007 

4 

1 

49.2212 

-68.117 

166.5589 

0.4100 

3 

1 

43.9908 

-76.742 

164.7233 

0.4741 

5 

6 

42.3197 

-73.303 

157.9424 

0.4721 

5 

7 

23.7103 

-93.996 

141.4164 

0.6923 

7 

6 

18.6095 

-97.013 

134.2321 

0.7518 

1 

2 

6.8621 

-114.990 

128.7144 

0.9119 

4 

3 

5.2304 

-108.991 

119.4518 

0.9283 

Comparisons  for  all  pairs  using  Tukey-Kramer  HSD 

q*  Alpha 

2.96474  0.05 


Abs{Dif)-LSD 

5 

7 

6 

4 

3 

1 

2 

5 

-177.46 

-153.75 

-132.00 

-58.90 

-58.92 

-19.49 

-9.80 

7 

-153.75 

-177.46 

-155.71 

-82.61 

-82.63 

-43.20 

-33.51 

6 

-132.00 

-155.71 

-171.12 

-97.98 

-98.10 

-58.75 

-49.01 

4 

-58.90 

-82.61 

-97.98 

-166.79 

-166.98 

-127.69 

-117.91 

3 

-58.92 

-82.63 

-98.10 

-166.98 

-177.46 

-138.03 

-128.34 

1 

-19.49 

-43.20 

-58.75 

-127.69 

-138.03 

-186.48 

-176.85 

2 

-9.80 

-33.51 

-49.01 

-117.91 

-128.34 

-176.85 

-180.91 

Positive  values  show  pairs  of  means  that  are  significantly  different. 


Level 

Mean 

5 

A 

162.6380 

7 

A 

138.9278 

6 

A 

120.3183 

4 

A 

49.3245 

3 

A 

44.0941 

1 

A 

0.1033 

2 

A 

-6.7588 

Levels  not  connected  by  same  letter  are  : 

Level 

-  Level 

Difference 

Lower  CL 

5 

2 

169.3968 

-9.797 

5 

1 

162.5347 

-19.490 

7 

2 

145.6865 

-33.507 

7 

1 

138.8244 

-43.201 

6 

2 

127.0771 

-49.006 

6 

1 

120.2150 

-58.749 

5 

3 

118.5439 

-58.918 

5 

4 

113.3135 

-58.895 

7 

3 

94.8336 

-82.629 

7 

4 

89.6032 

-82.605 

6 

3 

76.2242 

-98.097 

6 

4 

70.9938 

-97.976 

4 

2 

56.0833 

-117.909 

3 

2 

50.8529 

-128.341 

4 

1 

49.2212 

-127.686 

3 

1 

43.9908 

-138.034 

5 

6 

42.3197 

-132.001 

5 

7 

23.7103 

-153.752 

7 

6 

18.6095 

-155.712 

1 

2 

6.8621 

-176.851 

4 

3 

5.2304 

-166.978 

Upper  CL 
348.5905 
344.5598 
324.8802 
320.8495 
303.1605 
299.1790 
296.0062 
285.5221 
272.2960 
261.8118 
250.5453 
239.9635 
230.0755 
230.0466 
226.1281 
226.0159 
216.6409 
201.1726 
192.9306 
190.5756 
177.4390 


different. 


Difference 


■  I 


I 

I 
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Tests  that  the  Variances  are  Equal 


Level 

Count 

Std  Dev 

MeanAbsDif  to  Mean 

MeanAbsDif  to  Median 

1 

48 

286.1385 

216.4796 

213.5568 

2 

51 

274.5145 

192.9826 

192.0040 

3 

53 

294.5859 

238.2328 

236.6768 

4 

60 

278.4294 

228.2225 

226.7965 

5 

53 

345.1214 

279.0342 

278.9815 

6 

57 

308.7574 

247.6645 

247.8394 

7 

53 

359.1178 

297.0261 

294.7004 

Test 

F  Ratio 

DFNum 

DFDen 

Prob  >  F 

0'Brien[.5] 

1.2329 

6 

368 

0.2885 

Brown-Forsythe 

1.9188 

6 

368 

0.0768 

Levene 

1 .9868 

6 

368 

0.0667 

Bartlett 

1.1944 

6 

0.3057 

Welch  Anova  testing  Means  Equal,  allowing  Std  Devs  Not  Equal 

F  Ratio  DFNum  DFDen  Prob  >  F 
2.4046  6  162.36  0.0297 
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Oneway  Analysis  of  45km  E/W  By  45km  LT  Bins  (Nominal  CARP) 


Means  and  Std  Deviations 


Level 

Number 

Mean 

Std  Dev 

Std  Err  Mean 

Lower  95% 

Upper  95% 

1 

48 

-15.632 

397.339 

57.351 

-131.0 

99.74 

2 

51 

-10.805 

359.155 

50.292 

-111.8 

90.21 

3 

53 

11.717 

378.416 

51.979 

-92.6 

116.02 

4 

60 

5.949 

363.862 

46.974 

-88.0 

99.94 

5 

53 

65.746 

413.535 

56.803 

-48.2 

179.73 

6 

57 

-12.650 

442.575 

58.620 

-130.1 

104.78 

7 

53 

10.295 

441.901 

60.700 

-111.5 

132.10 

Means  Comparisons 

Comparisons  for  each  pair  using  i 

Student's  t 

t 

Alpha 

1 .96643 

0.05 

Abs{Dif)-LSD 

5 

3  7 

4  2 

6 

1 

5 

-153.14 

-99.11  -97.69 

-88.81  -78.08 

-72.03 

-75.70 

3 

-99.11 

-153.14  -151.72 

-142.84  -132.11 

-126.06 

-129.73 

7 

-97.69 

-151.72  -153.14 

-144.26  -133.53 

-127.48 

-131.15 

4 

-88.81 

-142.84  -144.26 

-143.93  -133.39 

-127.21 

-131.08 

2 

-78.08 

-132.11  -133.53 

-133.39  -156.11 

-150.10 

-153.71 

6 

-72.03 

-126.06  -127.48 

-127.21  -150.10 

-147.67 

-151.45 

1 

-75.70 

-129.73  -131.15 

-131.08  -153.71 

-151.45 

-160.92 

Positive  values  show  pairs  of  means  that  are  significantly  different. 


Level  Mean 

5  A  65.74553 

3  A  11.71682 

7  A  10.29480 

4  A  5.94860 

2  A  -10.80468 

6  A  -12.65008 

1  A  -15.63213 


Levels  not  connected  by  same  letter  are  significantly  different. 


Level 

-  Level 

Difference 

Lower  CL 

Upper  CL 

p-Value  Difference 

5 

1 

81.37766 

-75.698 

238.4538 

0.3090  1  1 

II 

5 

6 

78.39561 

-72.032 

228.8237 

0.3061  1  1 

1 1 

5 

2 

76.55021 

-78.083 

231.1830 

0.3310  1  1 

II 

5 

4 

59.79693 

-88.808 

208.4020 

0.4293  1  1 

1  1 

5 

7 

55.45073 

-97.688 

208.5895 

0.4769  1  1 

1  1 

5 

3 

54.02871 

-99.110 

207.1675 

0.4883  1  1 

1  1 

3 

1 

27.34895 

-129.727 

184.4251 

0.7323  1  1 

■ 

1  1 
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Levei 

-  Levei 

Difference 

Lower  CL 

Upper  CL 

p-Vaiue  Difference 

7 

1 

25.92693 

-131.149 

183.0030 

0.7457  n - 

■ 

- 1 - 1 

3 

6 

24.36690 

-126.061 

174.7950 

0.7503  1  1 

■ 

1  1 

7 

6 

22.94488 

-127.483 

173.3730 

0.7644  1  1 

■ 

1  1 

3 

2 

22.52150 

-132.111 

177.1543 

0.7747  1  1 

■ 

1  1 

4 

1 

21.58074 

-131.079 

174.2402 

0.7812  1  1 

■ 

1  1 

7 

2 

21.09948 

-133.533 

175.7323 

0.7886  1  1 

■ 

1  1 

4 

6 

18.59869 

-127.212 

164.4089 

0.8021  1  1 

■ 

1  1 

4 

2 

16.75328 

-133.391 

166.8976 

0.8264  1  1 

■ 

1  1 

3 

4 

5.76822 

-142.837 

154.3733 

0.9392  1  1 

l 

1  1 

2 

1 

4.82746 

-153.706 

163.3605 

0.9523  1 

l 

1  1 

7 

4 

4.34620 

-144.259 

152.9513 

0.9542  1 1 

1 

1  1 

6 

1 

2.98205 

-151.453 

157.4167 

0.9697  II 

1 

1  1 

2 

6 

1.84541 

-150.103 

153.7942 

0.9810  II 

1  1 

3 

7 

1 .42202 

-151.717 

154.5608 

0.9854  II 

1  1 

Comparisons  for  all  pairs  using 

Tukey-Kramer  HSD 

q* 

Alpha 

2.96474 

0.05 

Abs{Dif)-LSD  5 

3 

7 

4 

2 

6 

1 

5 

-230.88 

-176.85 

-175.43 

-164.25 

-156.59 

-148.40 

-155.44 

3 

-176.85 

-230.88 

-229.46 

-218.28 

-210.61 

-202.43 

-209.47 

7 

-175.43 

-229.46 

-230.88 

-219.70 

-212.04 

-203.85 

-210.89 

4 

-164.25 

-218.28 

-219.70 

-217.00 

-209.62 

-201.24 

-208.58 

2 

-156.59 

-210.61 

-212.04 

-209.62 

-235.37 

-227.24 

-234.19 

6 

-148.40 

-202.43 

-203.85 

-201.24 

-227.24 

-222.63 

-229.85 

1 

-155.44 

-209.47 

-210.89 

-208.58 

-234.19 

-229.85 

-242.61 

Positive  vaiues  show  pairs  of  means  that  are  significantiy  different. 


Levei  Mean 

5  A  65.74553 

3  A  11.71682 

7  A  10.29480 

4  A  5.94860 

2  A  -10.80468 

6  A  -12.65008 

1  A  -15.63213 


Leveis  not  connected  by  same  ietter  are  significantiy  different. 


Level 

-  Level 

Difference 

Lower  CL 

Upper  CL  Difference 

5 

1 

81.37766 

-155.442 

318.1971  1 - 1 - 

- n 

5 

6 

78.39561 

-148.401 

305.1920  1  1 

— 

1 1 

5 

2 

76.55021 

-156.585 

309.6859  1  1 

II 

5 

4 

59.79693 

-164.251 

283.8448  1  1 

1  1 

5 

7 

55.45073 

-175.432 

286.3339  1  1 

— — 

1  1 

5 

3 

54.02871 

-176.854 

284.9118  1  1 

1  1 

3 

1 

27.34895 

-209.470 

264.1683  1  1 

■ 

1  1 

7 

1 

25.92693 

-210.892 

262.7463  1  1 

■ 

1  1 

3 

6 

24.36690 

-202.429 

251.1633  1  1 

■ 

1  1 

7 

6 

22.94488 

-203.851 

249.7413  1  1 

■ 

1  1 

3 

2 

22.52150 

-210.614 

255.6572  1  1 

■ 

1  1 

4 

1 

21.58074 

-208.580 

251.7412  1  1 

■ 

1  1 

7 

2 

21.09948 

-212.036 

254.2352  1  1 

■ 

1  1 

4 

6 

18.59869 

-201.235 

238.4328  1  1 

■ 

1  1 

4 

2 

16.75328 

-209.615 

243.1217  1  1 

■ 

1  1 

3 

4 

5.76822 

-218.280 

229.8161  1  1 

j 

1  1 

2 

1 

4.82746 

-234.189 

243.8435  1 

j 

1  1 

7 

4 

4.34620 

-219.702 

228.3941  1  1 

j 

1  1 

6 

1 

2.98205 

-229.855 

235.8189  II 

j 

1  1 

2 

6 

1.84541 

-227.244 

230.9345  II 

j 

1  1 

3 

7 

1 .42202 

-229.461 

232.3051  II 

j 

1  1 

100 


Tests  that  the  Variances  are  Equal 


Level 

Count 

Std  Dev 

MeanAbsDif  to  Mean 

MeanAbsDif  to  Median 

1 

48 

397.3390 

326.6421 

326.6421 

2 

51 

359.1552 

299.1003 

298.4077 

3 

53 

378.4157 

295.0329 

294.6817 

4 

60 

363.8620 

294.2648 

294.1885 

5 

53 

413.5354 

338.0630 

335.4860 

6 

57 

442.5745 

366.6649 

360.2060 

7 

53 

441.9014 

362.9761 

362.5760 

Test 

F  Ratio 

DFNum 

DFDen 

Prob  >  F 

0'Brien[.5] 

1 .0776 

6 

368 

0.3753 

Brown-Forsythe 

0.8953 

6 

368 

0.4983 

Levene 

1 .0443 

6 

368 

0.3961 

Bartlett 

0.7998 

6 

0.5699 

Welch  Anova  testing  Means  Equal,  allowing  Std  Devs  Not  Equal 

F  Ratio  DFNum  DFDen  Prob  >  F 
0.2408  6  162.4  0.9624 


101 


Bivariate  Fit  of  45km  OP  N/S  By  45km  OP  E/W  (Nominai  OP) 


Linear  Fit 

Bivariate  Normal  Eiiipse  P=0.95(l 


Linear  Fit 

45km  OP  N/S  =  -1 .500542  +  0.0839099  45km  OP  E/W 

Summary  of  Fit 

RSquare 
RSquare  Adj 
Root  Mean  Square  Error 
Mean  of  Response 
Observations  (or  Sum  Wgts) 


Lack  Of  Fit 

Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Lack  Of  Fit 

371 

1250996.5 

3371.96 

Pure  Error 

2 

0.0 

0.00 

Prob  >  F 

Totai  Error 

373 

1250996.5 

Max  RSq 
1.0000 

Anaiysis  of  Variance 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Modei 

1 

9088.7 

9088.70 

2.7099 

Error 

373 

1250996.5 

3353.88 

Prob  >  F 

C.  Totai 

374 

1260085.2 

0.1006 

Parameter  Estimates 

Term  Estimate 

intercept  -1.500542 

45km  OP  E/W  0.0839099 

Correiation 

Variabie  Mean 


Std  Error 

t  Ratio 

Prob>|t| 

3.102866 

-0.48 

0.6290 

0.050973 

1.65 

0.1006 

Std  Dev 

Correiation 

Signif.  Prob 

0.007213 

0.004551 

57.91268 

-0.13898 

375 


Number 
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Variable 
45km  OP  E/W 
45km  OP  N/S 


Mean 

16.2265 

-0.13898 


Std  Dev 
58.74914 
58.04491 


Correlation 

0.084928 


Signif.  Prob 
0.1006 


Number 

375 
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Appendix  E:  Bivariate  Normal  MATLAB  Code 


function  [M]  =  BivarStats ( Z ) 

D  =  []; 

k  =  length ( Z ) ; 

X  =  Z  (  :  ,  2  )  ; 
y  =  Z  (  :  ,  1 )  ; 

Zbar  =  [mean(x)  mean(y)]'; 

Sigma  =  cov([x,y]) 

V  =  [  ]  ; 

for  i  =  1 :  k 

xt  =  (i  -  . 5) /k; 

Zi  =  Z (i,  :  )  '  ; 

di2  =  (Zi  -  Zbar) ' * (Sigma^-1) * (Zi  -  Zbar); 
Xsq  =  chi2inv (xt, 2 ) ; 

D  =  [D; di2 , Xsq] ; 

end 

M  =  sort (D) ; 
plot(M(:,l) ,M(:,2)  ,  '*') 


mu  =  [mean(x)  mean(y)] 
xl  =  -1250:10:1250; 
x2  =  -1250:10:1250; 

[XI, X2]  =  meshgrid (xl , x2 ) ; 

F  =  mvnpdf ( [Xl ( : )  X2 ( : ) ] ,mu, Sigma) ; 

F  =  reshape (F, length (x2 ),  length (xl ))  ; 
h  =  surf (xl , x2 , F) ; 

caxis ( [min (F ( : ) ) - . 5* range (F ( : ) ) , max (F (:))]) ; 
axis([-1250  1250  -1250  1250  0  2.0e-6]) 

xlabel  (  '  xl ' )  ;  ylabel  (  '  x2  '  )  ;  zlabel  (' Probability  Density'),' 
[V,D]  =  eigs (Sigma) 


Note:  This  script  is  optimized  to  view  CARP  data.  In  order  to  view  OP  data,  variables 
xl,  x2,  and  axis  must  be  changed  to  smaller  scales  to  ensure  correct  plotting  of  data. 
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Appendix  F:  Bivariate  Normal  MATLAB  Input  Script 


ooooooooooooooooooooooooooooooooooooo 

g, 
o 

Bivariate  Normal  Test  Script  % 

q, 
o 

ooooooooooooooooooooooooooooooooooooo 


r  =  menu  ('Choose  resolution  to  be  evaluated' 5km' ISkm' 45k;m' Skm 
OP' ,  ' 15km  OP' ,  ' 45km  OP' )  ; 

format  compact 

if  r  ==  1 


[150.437523 

LO 

1 

750732 

-131 . 954659 

-92  . 

498073 

-29.42237 

946. 

396638 

159.860508 

5.546833 

979.366917 

469. 

531369 

32.588577 

325. 

299697 

-46.562242 

403. 

046916 

32.588205 

-303 

.228644 

-46.558575 

73.977157 

-263.907812 

-90. 

609921 

-259.151387 

146. 

072017 

-164.545424 

-999 

. 964325 

285.463218 

-210 

.725415 

-37.245641 

-358 

. 683511 

-392 . 646767 

44.372029 

-136.484016 

-249 

.547516 

-364.766941 

CN 

1 

471685 

-214 . 158804 

-378 

.866974 

142 .817998 

-107 

.249294 

173.769628 

452  . 

84694 

-339.868473 

-210 

.723391 

108.586868 

-565 

. 640328 

166.076721 

107  . 

251901 

648.585893 

-232 

. 890433 

-79.142105 

258  . 

753925 

-378.584711 

-328 

. 952473 

-130.384523 

175. 

57118 

-104 . 021534 

380. 

865298 

325.91052 

-208 

.  949223 

69.833144 

199. 

638626 

-153.625722 

1112 

.763332 

-51.213623 

517  . 

617579 

436.088357 

700. 

628178 

121 . 055473 

-364 

.117126 

62 . 00744 

-81 . 

407993 

-249.822042 

-304 

. 999875 

333.765964 

317  . 

98451 

68.341809 

-260 

. 63892 

375.553112 

194  . 

099854 

105 


229.691587  -192.204881 

-296.493919  14.755564 
73.008638  -40.592793 

339.818724  11.096987 

153.662914  -515.7301 

-231.196065  303.120295 
-231.196065  303.120295 
-341.325626  386.306502 
260.718193  534.257617 

259.144239  245.780246 

310.303691  183.007095 

-142.764426  -299.455867 
215.643433  229.25381 

333.656969  221.82576 

-141.274091  110.910966 
302.64853  641.508008 

83.812547  123.775894 

-277.804449  317.982626 
353.822272  534.149857 

-302.586035  -203.293986 
-162.938233  -314.318438 
-330.558394  -351.246922 
116.366124  391.958386 

-243.590615  449.188988 
299.552512  -218.155673 

311.940967  260.643756 

-82.219532  -105.364325 

43.4882  33.273105 
266.866183  597.032814 

-6.146016  844.690329 

200.209338  -147.841063 

198.683175  -131.204918 

-103.958339  231.027973 
819.384442  20.331135 

-18.624713  253.207526 

530.763107  328.973875 

228.128984  -343.818945 

625.533231  562.000646 

10.802642  251.321943 

419.04352  839.26493 

499.732322  -512.058407 

-321.227325  -268.064812 
169.167677  40.705575 

-206.357891  478.80067 
74.500297  449.185559 

550.821917  64.787228 

-159.888899  -369.661988 
69.830193  -260.638981 

296.392731  -79.407124 

509.042745  -388.171497 

443.903561  759.743693 

563.420716  3.676363 

-152.067713  177.457235 
245.296661  64.663595 

133.525494  149.729633 
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276.213843 

-55.866525 

-100.85067 

-496.570319 

-319.646532 

139.68508 

-52.789655 

144.324797 

254.570803 

-85.296202 

133.439541 

228 . 141097 

628 . 60283 

622.379744 

170.701857 

-229.661101 

-180.082911 

457 . 85284 

479.608254 

239.005926 

-100.8314 

38.725473 

-169.201363 

-119.469029 

6.237754 

681.414155 

374 . 031016 

68.354335 

437.59987 

389.485016 

-167 . 612399 

-82.229343 

333.630493 

426.693836 

532.2713 

6.238602 

-122.523509 

141.169977 

-336.694439 

313.465104 

-74.49009 

460.975672 

131 . 910026 

451.494029 

159.875421 

114 . 901424 

97.788801 

-66.769215 

232.779182 

417.426942 

-310.357214 

204 . 856675 

35.666135 

251.427883 

-26.4407 


-205 

. 179526 

-96. 

0481 

305. 

003077 

-497 

.197876 

-234 

. 680845 

582  . 

27833 

-423 

.344789 

395. 

506704 

-654 

.367091 

284  . 

706548 

615. 

551626 

-146 

.065926 

-463 

.  9162 

-495 

.304997 

207  . 

070349 

-171 

. 90861 

-404 

. 820562 

-554 

.540096 

-114 

. 66905 

-158 

.  931519 

-428 

.889908 

179. 

342452 

40.5 

9455 

94.2 

74332 

341 . 

936633 

-79. 

498541 

136. 

759471 

266. 

184125 

-11 . 

081169 

672  . 

899756 

329. 

071647 

-430 

. 663267 

-321 

. 633801 

46.2 

5899 

-926 

.085193 

-323 

.525046 

-600 

. 800893 

377  . 

09516 

-365 

.997915 

149. 

734112 

1.77 

4849 

-432 

.538056 

240. 

344032 

22 . 1 

92612 

-249 

.546644 

-654 

.369744 

134  . 

977785 

671 . 

005521 

377  . 

097403 

-55. 

446123 

72 . 0 

96539 

414  . 

140448 

545. 

344413 

-556 

.432674 

253. 

208408 
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-214 . 159621 

678.328171 

529.216524 

274.701132 

367.758551 

-375.64323 

442.351826 

392.503062 

524.588187 

49.722655 

49.722655 

-189.294085 

-235.866458 


-297 . 680654 

410.390456 

-12 . 851172 

299.461054 

37 . 051021 

83.18982 

-399.265895 

279.169439 

-373.308803 

94.273773 

94.273773 

-133.09044 

138 . 640492]  ; 


disp('5k;m  Resolution') 
elseif  r  ==  2 


[153.603343 

-118 

.340326 

-130.371625 

-96. 

158125 

-4 . 655455 

857  . 

668393 

96.268581 

158  . 

934998 

952 . 920921 

462  . 

097767 

41 . 899782 

282  . 

821114 

-17.041904 

332  . 

729886 

44 . 971688 

-295 

.797607 

-97.772951 

79.523043 

-299.572521 

LO 

[ 

335875 

-282.431785 

127  . 

550683 

-155.233142 

-990 

.758945 

302.59606 

-253 

.203541 

-21.788671 

-345 

.707068 

-423.652789 

38 . 827821 

-122.519939 

-317 

.  979425 

-283.994818 

199. 

642461 

-336.829662 

1 

o 

698133 

-206.338611 

-438 

. 093185 

119.542256 

-81 . 

40749 

176.935344 

480. 

574498 

-319.662536 

-208 

.  949514 

113.24295 

-539 

.798258 

180.041702 

46.2 

51358 

660.968649 

-251 

.300742 

-17 . 038888 

225. 

480539 

-373.93228 

-408 

.475292 

-66.682078 

218  . 

159926 

-94.708016 

476. 

913246 

335.222758 

-223 

. 699951 

88.456023 

122  . 

001625 

-116.384939 

961 . 

259364 

-74.49277 

488  . 

004722 

402.008993 

624  . 

875153 

110.160146 

-334 

. 615175 

145.802611 

-170 

.  135359 

-251.404531 

-288 

.363305 

108 


311 . 973734 

104 . 002275 

397.254163 

263.863261 

-318.191049 

114 . 914741 

304 . 067005 

153.662125 

-209.5029 

-243.565856 

239.021988 

268.456026 

328 . 926093 

-172.286067 

190.876823 

330.585192 

-159.899366 

240.628771 

110.168341 

-184 . 673758 

232.777765 

-266.834078 

-69.830564 

-338.379605 

246.697641 

-311 . 939721 

311 . 93883 

313.52413 

-27 . 934305 

32.592857 

268.448464 

-72 . 914363 

217.342994 

270.001954 

128.71596 

665.746179 

-63.603479 

554 . 041767 

239.020896 

637 . 918989 

91.542849 

439.244126 

450.093585 

-363.126678 

166.095351 

-207 . 940765 

80.739958 

543.093614 

-32.592013 

294.779829 

316.599272 

470.300658 

493.533006 

450.082195 

-166.035713 


336.394839 

-244.002064 

94.281662 

-262.410237 

9.210747 

-94.272833 

40.708876 

-465.820635 

218.051675 

336.39392 

563.759167 

242 . 120465 

114.575995 

-312.431852 

192.209274 

184 . 892615 

122 . 00225 

645.166305 

0.000625 

395.617355 

536.031675 

-194.089475 

-299.457576 

-338.270164 

334 . 620084 

380.759388 

-279.155861 

255.0983 

-157 . 159683 

36.933093 

621.100364 

863.212577 

-123.773211 

-181.112851 

140.414051 

73.99993 

240.342155 

336.406154 

-245.884993 

563.77601 

275.389831 

985.334384 

-495.424346 

-271.723373 

36.93457 

488.006237 

417.798085 

72 . 106872 

-260.638659 

-369.659564 

-81.291958 

-378 . 857016 

830.063075 

57.350868 

192.208499 


109 


246.787913 

45.067885 

181.595579 

-41 . 899711 

-125.713893 

-654 . 859256 

-372.440098 

94 . 613073 

-77 . 648291 

239.026172 

273.193393 

-82.22249 

307.300395 

265.387662 

656.543356 

355.417472 

189.326184 

-102.40256 

-246.753778 

515.210714 

465.639347 

220.384254 

-135.00049 

147.362972 

-148.994413 

86.969034 

96.268581 

731 . 05113 

369.375483 

155.240349 

381.738822 

498 . 157502 

-170.7763 

-71.333709 

44 . 974787 

397 .188083 

533.855149 

3.072753 

-13.965528 

-21 . 69566 

159.792038 

-363.140294 

356.945776 

237.473086 

-142.739924 

481.092259 

200.148247 

515.1713 

74.491026 

156.802296 

34 . 179215 

-83.810362 

367.793351 

397.221654 

-159.783525 


16.639625 

CO 

( — 1 

1 

85368 

-123 

.884751 

-53. 

569553 

323. 

525345 

-456 

.48445 

-240 

.224476 

611 . 

89076 

-425 

.230098 

164  . 

482681 

-654 

.366585 

380. 

865587 

373. 

327981 

-110 

. 90648 

-502 

.732839 

-441 

.748828 

262  . 

525673 

-168 

.25074 

-447 

.4086 

-521 

.264182 

-120 

.215227 

-147 

. 84093 

-428 

.889492 

88.72956 

16.637649 

338.276299 

158 . 934998 

-79.494933 

134 . 98473 

288.367119 

-81.400653 

449.198854 

449.187355 

-428.888793 

-389.960081 

-151.495382 

-940.947056 

-351.252589 

-658 . 142308 

57.229667 

495.436504 

-417 . 681113 

59.121806 

184.778867 

118.231327 

-458.490012 

120.229077 

-27.713962 

-306.888238 

-656.143721 

219.934316 

713.595102 

317 . 875527 

-62.767043 

336.391785 


no 


187 .81679 

390.072607 

27 . 937069 

493.549411 

310.373693 

-554.545499 

18.620573 

64 . 660754 

-172.258202 

-262.522956 

667.431061 

430.686227 

459.389208 

-64 . 649835 

263.811928 

88.731619 

347 . 654763 

-145.951772 

-327.498562 

142.303154 

274.722248 

-345.591659 

439.055548 

207.079655 

491 .99472 

-397.377965 

-26.351316 

31.387613 

-248.2361 

-247.770277 

-145.822316 

110.911195] ; 

disp('15k;in  Resolution') 

:if  r  ==  3 

Z  =  [41.89052 

183.002514 

-131 . 951595 

133.09296 

9.311074 

687 . 643011 

-152 . 047253 

-475.138497 

808 . 612362 

88.761618 

-9.311145 

144 .183257 

-52.792419 

-216.274911 

-378.704418 

151 . 621515 

-262.3174 

181.119679 

-189.314262 

20.298385 

-99.360253 

-977.783201 

89.942178 

-395.505738 

-27 . 933667 

-162.705245 

-628.51875 

-314.299165 

-49.622428 

-316.094594 

-142.748589 

-227.364732 

-136.504051 

-456.505509 

-41.895172 

29.502256 

209.531727 

340.052183 

-296.38337 

-192.313728 

319.688816 

-3.765681 

124 .184107 

190.434181 

606.715086 

-576.71563 

-186.224339 

-138 . 635981 

-88.455551 

-458.39169 

429.965892 

709.832879 

127.303763 

321 . 639652 

46.558893 

-260.638691 

-41 . 901027 

-64 . 660549 

-215.641832 

789.350159 

108 . 669223 

229.25174 

513.761038 

525.061397 

-173.834728 

-587 . 823039 

400.39536 

-482.451074 

Ill 


-439.213951  -308.764102 
9.310746  -181.116495 

356.94296  -182.994982 

49.719154  -338.276337 

-263.899927  88.731646 
321.263092  327.189956 

336.8337  741.216861 

111.719394  170.137199 

175.361712  -500.867803 

-161.476445  517.507249 
-265.379784  -127.54305 
-526.083054  -271.715804 
-32.593824  377.093893 

-198.61166  521.279658 

-45.065089  -301.343094 

-9.311655  225.480242 

190.852911  332.733068 

-209.479928  97.93611 
-121.06031  319.865163 

71.325144  -133.091962 

12.382893  345.707426 

48.146046  264.409444 

117.977183  225.481011 

304.089758  -341.931708 

259.256891  510.189294 

97.773676  -118.229693 

110.168698  597.028684 

-116.403835  33.273702 
437.635693  236.581216 

-200.175234  -345.594602 
-68.247284  -234.79689 

-555.626012  -500.965346 
403.369767  434.444804 

-99.264751  112.796059 

211.09252  -205.181512 

-86.888732  -68.431042 

1.582977  247.773172 

9.31141  -184.887259 
-144.338064  -613.775408 
-24.863944  -61.000478 

-136.598673  619.212142 
58.947082  145.957742 

-37.249681  611.779359 

-287.059573  -53.565401 
-111.747384  175.571187 
-861.317568  -774.447802 
-38.827538  -404.822315 

-68.353151  499.095061 

-329.020321  -208.949618 
-43.391964  79.411753 

-170.671229  -380.864347 
-21.693814  -428.890231 

-265.393168  46.25312 
484.233777  316.105606 

286.938478  27.731967 


112 


27 . 935719 

391 . 13757 

-41 . 900652 

130.372965 

274 . 635772 

-526.079449 

54.286875 

782 . 151003 

-166.002438 

-324.270791 

4 . 656204 

229.566216 

-180.100768 

-63.59445 

266.931805 

996.41853 

-4 . 655514 

-21 . 698092 

46.554623 

-85.388264 

-315.122406 

-293.327171 

-201.794848 

366.188484 

7.728465 

-487.255053 

-72 .898168 

-297 . 964104 

-332 . 175166 

-198 . 611054 

161.436953 

194 . 056139 

152.07211 

-277.760738 

135.035257 

304 . 133344 

724 . 654174 

-88.455214 

-290.13476 

-114 . 829983 

-186.250033 

359.994121 

20.112171 

-175.324809 

-130.362645 

346.139233 

268.422354 

352.249842 

-270.028221 

299.589733 

162 . 940434 

420.540572 

383.31989 

-400.418069 

43.488936 


-171 . 910615 

-635.839085 

-255.093407 

571.187194 

-547.228539 

-110.895974 

-567.527222 

-120.084244 

280.9373 

-138 . 632585 

146.068447 

-77 . 634662 

-186.770817 

768 . 940409 

-33.269429 

369.71408 

266.184491 

364.228332 

-312.434138 

92.3885 

-12 . 971365 

114.574527 

-1042.55308 

-271.723371 

-404 . 822209 

-116.443502 

-406.70824 

123.780315 

-90.607793 

-478 . 686474 

-24.066213 

244.004156 

427.004546 

-201.519822 

-574 . 956431 

85.073139 

-373.29774 

9.205952 

-476.909765 

305.002817 

464 . 049011 

-281.039591 

90.613654 

820.736947 

465.823284 

523.16851 

-181.223775 

-345.700709 

-329.066844 

277.279443 

-269.843345 

-48.014943 

-293.90478 

-40.695742 

-123.775476 


113 


327.519001 

380.127918 

175.361712 

-82.325084 

-226.545833 

-313.40554 

79.155737 

239.020624 

211 .19562  A 

-21 A . 673581 

226.554153 

490.267357 

133.540639 

10.893018 

-111.745279 

18 . 622338 

-100.858404 

384 . 878957 

150.462335 

-77.5746 

-40.32004 

366.191412 

27 . 936947 

231.208915 

308 .806255 

181.592677 

71.333994 

257.584089 

-296.429392 

-63.601676 

229.706804 

-229.595072 

27 . 932369 

-423.674402 

552.320702 

-387 . 92265 

229.719789 

-3.073298 

276.278688 

-4 . 655707 

-226.558421 

-71.42489 

296.380548 

386.46779 

-273.12814 

-99.348964 

-18 . 624203 

-788.391246 

-239.014115 

-167 . 619385 

-158.275651 

-299.455235 

-85.289644 

439.172981 

-257.565563 


805.987976 

101.712276 

-500.867803 

73.866363 

73.868815 

25.847146 

462 .162123 

619.21426 

-184 .882966 

-301.339314 

210.731833 

171 . 923476 

3.660948 


42.4 

78673 

554  . 

550683 

83.1 

82671 

99.8 

19448 

101 . 

712099 

-11 . 

089862 

164  . 

590666 

-35. 

158421 

-310 

.541955 

586. 

048323 

-11 . 

088261 

-96. 

043131 

306. 

889587 

96.048299 

240. 

345309 

347  . 

485214 

244  . 

002282 

463. 

939963 

-292 

.024261 

-349 

.367412 

-373 

.425337 

382  . 

65692 

380. 

762125 

-382 

. 637304 

-319 

.864197 

210. 

733121 

-188 

.547298 

-366 

. 000746 

40.7 

04242 

693. 

192788 

-1 . 8 

77784 

794  . 

784651 

-62  . 

774685 

364  . 

117627 

-565 

.  610155 

-267 

. 956236 

739. 

438279 

-410 

.367131 

247  . 

666978 

-447 

.411609 

-199 

. 628698 

253. 

32203 

114 


481 . 

050976 

355. 

260044 

206. 

347838 

350. 

809345 

-307 

.280802 

60.531261 

119. 

440033 

-127 

.288406 

195. 

549289 

333. 

762981 

634  . 

744934 

-37  . 

243053 

-183 

. 041729 

-114 

. 821726 

764  . 

819651 

-51 . 

2168 

256. 

053007 

591 . 

203069 

772  . 

884843 

248  . 

329163 

-17  . 

042449 

218  . 

804394 

3.07277  325 

-288 

. 671567 

-138 

. 101325 

-131 

. 948909 

518  . 

223072 

113. 

227609 

40.3 

18866 

128  . 

853161 

-76. 

076103 

-152 

.16538 

-204 

. 849697 

676. 

558157 

496. 

688858 

66.677766 

-245 

.16719 

99.2 

73893 

795. 

919582 

665. 

623738 

-367 

.797761 

-179 

. 985721 

1 

CO 

12798 

-221 

. 908641 

373. 

941534 

209. 

504494 

-32  . 

588564 

167  . 

614378 

257  . 

664028 

591 . 

206676 

291 . 

720483 

-147 

.399061 

392  . 

709391 

290. 

215881 

218  . 

829008 

652.495679 
-171 . 905346 
16.638709 
-530.476177 
-550.775155 
417.79798 
-282 . 821104 
-500.980506 
-805.982145 
92.393733 
0.020744 
-214.389531 
-412.252154 
221 . 820654 
286.512137 
22 .182144 
798.556609 
183.020069 
846.606965 
160.822959 
22 .182001 
-208 . 952512 
29942 

134 . 981835 

-219.933625 

-702.3931 

-110.896499 

-432.549092 

-447.411615 

-499.096006 

-149.728395 

-85.06678 

-480.571817 

62.798863 

205.196376 

323.413655 

-64.768456 

-545.343908 

349.400862 

-256.845949 

-262.517487 

-563.755126 

521.276468 

826.281871 

245.784135 

-35.04536 

728.347187 

160.821052 

828.168371 

7.338084 

-168.246415 

-388.184606 

460.284205 

-125.657024 

330.84719 
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521.390033 

-364.754657 

442.347029 

529.315646 

827 . 009115 

144.348939 

-26.447455 

-74.493181 

-754.203834 

246.771856 

257 . 647552 

370.935714 

-375.605624 

679.711464 

884.445878 

133.540639 

170.752051 

-296.410217 

231.202224 

-15.458803 

200.207624 

446.912466 

-76.084094 

-58 . 94495 

414.333381 

43.487364 

305.704068 

384 . 90422 

310.37511 

15.459112 

442.339117 

99.365395 

-277.785477 

282.506623 

-209.486469 

161.447066 

-324.41595 

802 . 17571 

-245.184592 

259.149663 

139.693558 

650.295392 

-32.590655 

-161.365601 

6.239364 

183.064379 

114 . 825002 

-308.788162 

-104 . 007145 

-263.816949 

-391 . 065752 

24.770797 

-32.586505 

-86.875589 

197 . 122178 


-164 

.466019 

CO 

1 

835731 

-293 

. 901427 

883. 

633874 

96.1 

94541 

-443 

. 638792 

395. 

504944 

-207 

. 068958 

-223 

.676646 

769. 

052964 

454  . 

734942 

-547 

.222619 

-330 

. 838059 

841 . 

055463 

271 . 

770752 

3.660948 

LO 

1 

793632 

362  . 

347467 

-164 

.47705 

454  . 

730785 

473. 

255184 

105. 

375056 

225. 

591133 

-255 

.093031 

-438 

.086834 

842  . 

916029 

-66. 

541262 

-447 

.292192 

578  . 

622391 

CO 

I 

04439 

354  . 

921987 

658  . 

029048 

55.4 

5898 

277  . 

279499 

-380 

.752978 

-227 

.364798 

-445 

.520754 

462  . 

197206 

12 . 979573 

-565 

. 638237 

-142 

.407346 

733. 

914187 

-345 

.706924 

208  . 

845179 

-266 

.184003 

530. 

484779 

-184 

.886251 

979. 

672763 

-310 

.548044 

175. 

574126 

-330 

. 837402 

591 . 

482888 

-512 

.073151 

242  . 

117241 

-756 

.073063 

116 


622  . 

308903 

-195 

.  958645 

-152 

.058568 

497  . 

322062 

-6.2 

35816 

-98  . 

045215 

274  . 

698099 

-60. 

996676 

-122 

.542636 

-770 

. 824941 

34 . 1 

77668 

194  . 

092446 

183. 

119278 

-427 

.003749 

-294 

. 904874 

-804 

.094242 

485. 

791246 

-624 

. 855759 

361 . 

712633 

-227 

.358643 

684  . 

382394 

-123 

.862572 

164  . 

519902 

-138 

. 636484 

-71 . 

318138 

-515 

.732963 

18 . 625205 

-3.660011 

662  . 

422451 

275. 

413548 

-347 

.720213 

-109 

.  12925 

180. 

079382 

567  . 

418011 

462  . 

443947 

138  . 

64894 

675. 

171913 

-40. 

680495 

316. 

581918 

105. 

369842 

45.072482 

353. 

137413 

384  . 

816848 

-220 

.038378 

195. 

537604 

491 . 

666647 

-223 

.479798 

476. 

915956 

-122 

.55434 

-580 

.391308 

198  . 

616157 

-428 

.776704 

555. 

46956 

-188 

.53166 

577  . 

344124 

-975 

.992152 

-124 

. 12601 

-720 

. 915286 

-66. 

662856 

-765 

.281477 

-243 

. 686327 

-177 

.453169 

-13. 

968208 

236. 

5711 

-335 

.205851 

-401 

. 045589 

572  . 

559931 

155. 

291457 

771 . 

337249 

813. 

33472 

89.96106 

112  . 

684913 

-460 

. 904623 

61.0 

11547 

180. 

013636 

-449 

.183818 

727  . 

686684 

321 . 

667862 

642  . 

353484 

-306 

.757242 

-363 

. 151632 

-517 

.  611034 

124  . 

213337 

-702 

.39374 

-121 
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613. 

125197 

646. 

957138 

474  . 

850912 

-186 

. 650402 

181 . 

589982 

-121 

.999335 

117 


504.352904 
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19.552358 
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30.911694 
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-4 .189841 
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-18.342093 
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20.204079 
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20.110984 
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37.242597 
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-43.295103 
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-6.238093 

-55.119935 
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13.779812 
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10.425591 
24.733005 
82.073614 
-92.386312 
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-65.658725 
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-26.618471 
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-4.214525 
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-27.061988 
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17 .191167 
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-47 .802325 
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i('5k;m  Resolution  OP') 
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-54.468002  14.085758 

-10.707359  81.075417 
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-6.610645  3.881862 

20.949185  -2.661824 
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24 . 02157 

-1.219983 

-10.800447 

-44.031378 

-30.259937 

95.604703 

-28.863301 

83.51548 

43.015169 
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210.175198] 

disp('45k;m  Resolution  OP') 

end 

[M]  =  BivarStats  f(Z); 
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Appendix  G:  JPADS-MP  Help 


JPADS-MP  Operation  -  N  Mission 

The  JPADS-MP  is  developed  by  Draper  Labs  and  Planning  Systems  Ine  (PSI), 
and  a  complete  user’s  manual  is  available  from  them.  The  following  sequence  of  figures 
will  provide  the  reader  with  sufficient  familiarity  with  the  JPADS-MP  Graphical  User 
Interface  (GUI)  to  recreate  the  steps  taken  in  this  research.  Upon  starting  the  JPADS- 


MP,  the  user  is  presented  with  the  main  GUI  page  as  shown  in  following  figure; 


Figure  45.  JPADS-MP  main  GUI. 
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In  these  examples,  the  N  mission  data  already  been  entered.  To  begin  anew,  the  user 
would  need  to  enter  a  mission  name  and  save  the  file.  Clicking  on  the  Aircraft  tab  allows 
the  user  to  selects  either  the  C-130  (the  default)  or  the  C-17.  The  next  figure  shows  the 
Drop  Zone  tab: 


Figure  46,  Drop  Zone  GUI 

The  Drop  Zone  tab  allows  the  user  to  select  the  desired  PI  Wind  Forecast  Reference  Point 
in  terms  of  Latitude,  Longitude,  and  Elevation.  This  is  the  point  that  the  wind  file  will  be 
centered  on.  For  this  study,  the  coordinates  for  Site  16  at  YPG  are  used  as  this  is  the 
point  from  which  the  weather  balloons  were  launched.  This  window  is  also  where  the 
aircraft  approach  data  is  defined.  It  is  worth  noting  that  the  window  is  titled  Release  Pass 
1.  A  key  feature  of  JPADS  is  the  ability  to  drop  on  multiple  Pis  or  make  multiple  release 
passes.  For  this  research,  only  a  single  Release  Pass  is  used. 

The  next  tab  is  Drop  Parameters. 
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[Drop  Parameters  -  N  -  Release  Pass  1 


3 
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|14:14 
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Figure  47,  Drop  Parameters  GUI 

This  window  allows  the  user  to  define  details  of  airdrop  in  terms  of  Altitude,  Airspeed, 
Date,  and  Time.  The  user  has  the  option  to  manually  specify  a  release  point,  but  that 
option  defeats  the  purpose  of  this  research  and  is  not  used.  It  is  critical  for  the  data  and 
time  set  here  be  in  agreement  with  the  date  and  time  of  the  weather  data  the  user  intends 
to  use  for  mission  planning,  otherwise  the  results  will  be  invalid.  The  user  must  next  fill 
out  the  Load  <&  Chute  tab.  For  full  details  of  options  for  this  page,  please  refer  to  the 
most  current  PSI  published  JPADS-MP  user’s  manual. 
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Figure  48,  Load  &  Chute  GUI 

For  this  research,  the  settings  in  Load  &  Chute  tab  are  held  constant.  The  guided 
parachute  system  is  the  Screamer  with  850  ft^  parachute.  The  Ballistic  Chute  Type  is  set 
to  two  G1 1  parachutes.  The  PI  coordinates  and  elevation  are  set  here.  This 
representative  mission  was  created  using  an  actual  test  point  from  the  ACTD  program, 
only  the  Run-In  heading  was  changed  to  a  cardinal  direction.  As  a  result,  the  PI  is  set  as 
being  the  JPADS  Center  PI  target  as  used  in  testing.  This  is  located  3.7  km  from  Site  16. 
It  may  have  been  better  for  the  purposes  of  this  analysis  to  have  set  Site  16  to  be  the  PI. 
Unfortunately,  this  was  realized  too  late  for  implementation.  However,  any  error 
incurred  by  this  is  believed  to  be  minimal  when  considering  that  the  minimum  tested 


135 


weather  resolution  was  5  km.  Total  Rigged  Weight  is  8000  lbs  at  Flight  Station  677  on 
the  left  side  of  the  fuselage.  CNUI  ID  is  set  as  1001,  but  is  not  needed.  The  Glide  Safety 
Factor  is  left  at  the  default  setting  of  0.89.  Finally  the  user  is  ready  to  gather  weather  data 
with  the  Weather  tab. 


Figure  49,  Weather  GUI 

The  Weather  tab  is  where  most  of  the  work  in  this  research  was  done.  The  Drop 
Information  shows  data  copied  from  the  previous  tabs.  The  next  step  is  to  acquire 
weather  data.  This  is  done  by  selecting  one  of  the  options  under  the  Weather  Acquisition 
section.  The  options  relevant  to  this  research  are:  Dropsondes,  4D  Forecast,  Balloon,  and 
Climatology.  In  the  interest  of  space,  only  the  4D  Forecast  method  will  be  shown  here. 
Other  methods  work  similarly;  where  there  are  differences,  they  are  detailed  in  the 
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manual.  Selecting  an  option  under  Weather  Acquisition  brings  up  the  Weather  Source  - 
Acquire  4D  Forecast  GUI. 


Figure  50.  Weather  Source  GUI. 

The  JPADS-MP  uses  the  4D  Forecasts  generated  by  AFWA.  These  come  in  a  format 
known  as  GRidded  Information  in  Binary  format  (GRIB)  files.  Once  these  are 
downloaded,  the  Browse  button  is  used  to  point  the  Weather  Source  GUI  to  the  location 
of  required  GRIB  files.  Once  the  appropriate  path  is  specified  in  the  “GRIB  Files 
Location”  field,  select  the  “Acquire  Forecast”  button.  This  will  read  the  weather  forecast 
into  the  JPADS-MP  Environmental  Data  folder.  Once  this  is  complete,  a  windgui 
window  opens  showing  the  duration  of  valid  times  for  the  forecast  as  well  as  the  forecast 
coverage  area. 
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Figure  51,  “windgui”  Information  Window 

The  user  ean  then  elose  this  window  as  well  as  the  Weather  Souree  window  and  return  to 

the  Weather  tab. 


Figure  52,  Weather  GUI  with  4D  Forecast  loaded. 

Here,  the  4D  Forecast  inventory  now  shows  an  increment  of  one  and  the  Wind  File 
Production  section  now  has  the  options  for  wind  file  generation  via  LAPS  Forecast-only, 
Forecast  Only,  and  Climatology  only.  These  are  listed  in  order  of  best  to  worst  methods 
for  calculating  an  accurate  wind  file.  The  Local  Analysis  and  Prediction  System  (LAPS) 
is  the  most  advanced  modeling  method  included  within  the  JPADS-MP.  It  allows 
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complex  modeling  of  wind  interaction  with  terrain  features  such  as  how  wind  will  flow 
over  or  around  terrain  obstacles.  Select  either  the  Best  Available  or  LAPS  Forecast-only 
buttons  to  begin  Wind  File  produetion. 


Figure  53,  Weather  GUI  during  Wind  File  production  and  production  complete. 

The  LAPS  Wind  File  produetion  takes  a  few  minutes  to  run.  During  this  time,  the  Build 
Status  indicator  will  be  yellow  and  display  “In  Progress”.  This  will  change  to 
“Complete”  and  the  status  bar  will  be  full  once  the  LAPS  Wind  File  has  been  generated. 
The  “Weather”  tab  can  now  be  closed  to  return  to  main  JPADS-MP  GUI  page.  Once 
here,  the  user  selects  the  Calculate  button  under  the  CARP/LAR  seetion.  Unless 
FalconView  is  installed,  do  not  select  generate  footprints  as  it  will  result  in  a  crash  of  the 
JPADS-MP  (this  option  was  not  used  in  this  research). 
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Figure  54,  JPADS  Main  GUI  CARP  Solution  TAB  after  successful  CARP 

calculation. 


Selecting  Compute  CARP  will  automatically  open  the  CARP  solution  tab.  The  CARP 
section  shows  the  Latitude  and  Longitude  of  the  Early,  Nominal,  and  Late  CARPs  which 
also  define  the  boundaries  of  the  LAR.  In  order  to  collect  this  data,  an  Optical  Screen 
Reader  tool  was  developed  by  Captain  Ryan  Eggert  of  the  Air  Eorce  Research  Eaboratory 
Advanced  Architecture  and  Integration  Branch. 
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Figure  55,  Screen  OCR  Main  GUI 

The  Screen  OCR  tool  reads  the  values  in  the  Early,  Nominal,  and  Late  CARP  coordinate 
boxes  and  copies  them  to  a  text  file.  In  doing  so,  it  also  converts  them  from  a  DDD 
MM.mmm  format  to  a  DDD.dddddd  format.  The  conversion  to  decimal  degrees  allows 
for  easier  mathematical  operations  later.  Additionally,  the  Screen  OCR  copies  the 
coordinates  for  the  Screamer  OP  from  its  memory  location  and  writes  it  to  the  same  text 
file.  To  use  the  Screen  OCR,  ensure  that  the  Mission  Name  is  correct.  This  must  match 
name  of  the  Screamer  mission  data  in  memory.  It  will  be  mission  followed  by 
“_R1  .scm”.  Check  the  type  of  data  from  which  the  wind  profile  is  being  generated  and 
type  in  the  launch  time  for  a  weather  balloon  (i.e.,  1414Z)  or  the  initialization  time  for  a 
weather  forecast  (i.e.,  0600  ini).  Select  “Specify  File.” 


Save  As 


Savein;|Q  20050G20  O  E  Cf HI’ 


Filename:  |20050620  1414Z  DalaM 

Save  as  type:  | 


,  I  ^^Cancel^^ 


Figure  56,  Screen  OCR  File  Save  As  GUI 
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If  starting  a  new  file  (for  a  weather  balloon),  enter  the  file  name  “YYYYMMDD  TTTTZ 
Data.txt”.  Otherwise,  click  the  appropriate  existing  file  and  select  Save. 

The  method  of  building  text  files  for  analysis  is  to  segregate  the  data  by  weather 
balloons.  The  Screen  OCR  allows  for  a  new  file  to  be  opened  and  then  to  append 
subsequent  data  to  this  file.  First,  the  CARP/LAR/OP  is  calculated  for  a  weather  balloon. 
This  data  is  saved  to  a  new  file  bearing  the  date  and  time  of  the  balloon  launch  as  the  file 
name.  Next  the  CARP/LAR/OP  is  calculated  for  each  weather  forecast  that  was  valid  for 
the  time  of  that  weather  balloon  launch.  Each  new  data  set  is  appended  to  the  text  file 
resulting  in  a  file  similar  to  the  one  shown  below: 


Figure  57,  Sample  text  file  record  of  CARP  and  OP  calculations  from  the  JPADS- 
MP  captured  by  the  Screen  OCR  program. 

As  can  be  seen,  each  line  represents  a  different  weather  input:  weather  balloon  on  the  first 

line,  followed  by  weather  forecasts  of  varying  resolution  and  initialization  time.  The 

coordinates  of  the  CARPs  and  OP  are  to  the  left  of  the  metadata.  Capt  Eggert  also 

developed  a  CARP  Analysis  tool  to  generate  Northing  and  Easting  data  from  the  raw 

coordinates  captured  by  the  Screen  OCR.  Note  that  in  the  Screen  OCR  created  text  file 

the  column  contain  information  in  the  following  order: 
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Early  Lat/Early  Lon  -  Nom  Lat/Nom  Lon  -  Late  Lat/Late  Lon  -  OP  Lat/OP  Lon  -  Metadata 
However,  the  CARP  Analysis  tool  changes  this  order  to: 

Nom  NS/Nom  EW  -  Early  NS/Early  EW  -  Late  NS/Late  EW  -  OP  NS/OP  EW  -  Metadata 


Figure  58,  CARP  Analysis  Tool  Initialization  GUI, 

Selecting  the  CARP  Analysis  tool  opens  the  GUI  above.  Selecting  Generate  Northing 
Easting  brings  up  the  Open  window  for  fde  selection. 


Figure  59,  File  selection  GUI  for  CARP  Analysis  Tool, 

After  finding  the  correct  file  to  analyze,  select  Open.  This  opens  the  Save  window. 
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Figure  60.  File  Save  As  GUI  for  CARP  Analysis  Tool. 

Be  extremely  eareful  to  change  the  file  name  here,  otherwise  you  will  destroy  you  data  as 

well  a  getting  no  results!  Northing/Easting  data  in  this  study  is  noted  by  an  “NE”  in  the 

file  name  to  distinguish  it  from  the  raw  data  text  file.  The  CARP  analysis  tool  functions 

by  comparing  each  weather  forecast  to  the  weather  balloon  data  in  line  one  of  the  text 


file.  This  results  in  a  file  similar  to  this: 


1^  20050620  1414Z  NE  Data.t><t  -  Notepad 


File  Edit  Format  View  Help 


jalxi 


Early  EW 

Late  NS  Late 

EW  OP  NS  OP  EW 

Source  Lead 

-Time 

67 

953. 587374 

462.204721 

950.764528 

462.101579 

79.047517 

80.632102 

15km  Data 

814 

33 

731.490987 

-81.268824 

729.028579 

-81.381192 

78.116919 

22.404192 

15km  Data 

2014 

69 

979.956292 

467.752795 

977.194776 

467.649755 

82. 585575 

79.855759 

5km  Data 

814 

41 

681.824147 

-81.272429 

679.421502 

-81.384795 

69.085537 

20. 518648 

5km  Data 

2014 

65 

774.933989 

844.825326 

772.419161 

844.728703 

98.878861 

116.456310 

45km  Data 

1414 

80 

998,602194 

367.936515 

995.818261 

307.831778 

86.961669 

70.650243 

45km  Data 

814 

8 

809.192097 

86.986301 

806. 543093 

86.876783 

66.292295 

24.400489 

45km  Data 

214 

95 

675. 578024 

-40. 569259 

673.182890 

-40.680933 

77.278890 

30. 500636 

45km  Data 

2014 

Figure  61,  Sample  text  file  containing  output  from  the  CARP  Analysis  Tool, 

In  this  file,  the  data  represents  error  in  the  forecasting.  A  value  in  the  Nominal  NS 
column  of  -40.578024  means  that  particular  forecast  generated  a  Nominal  CARP 
coordinate  that  was  40.578024  m  South  of  the  correct  Nominal  CARP  coordinate  as 
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defined  by  the  Nominal  CARP  calculated  from  the  weather  balloon  (an  actual  sampling 
of  the  atmosphere).  Also  note  that,  while  the  resolution  data  is  unchanged,  the 
initialization  time  has  been  replaced  by  the  Lead-Time.  This  is  accomplished  by  simply 
taking  the  difference  between  the  weather  balloon  launch  time  and  the  forecast 
initialization  time.  The  data  from  each  weather  balloon  (and  its  corresponding  forecasts) 
is  saved  in  a  folder  named  for  the  day  the  balloons  were  launched  on. 
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Appendix  H:  Weather  Balloon  Data  Listing 


Drop  Date  |  #  WxBs  |  WxB  File  |  WxB  Launch  Time 


20050620 


20050621 


20050624 


20050815 


20050816 


20050817 


20050818  3 


20050819  3 


20050912 


20050913 


20050915 


20051019 


20051020 


20051021 


20060125 


I  20060126  I  2 


1234Z 


1654Z 


146 


20060227 

2 

20060228 

2 

20060301 

3 

20060303 

2 

20060327 

3 

20060328 

3 

20060329 

1 

20060330 

3 

20060508 

2 

20060509 

2 

20060510 

2 

20060511 

2 

20060613 

2 

20060614 

2 

20060615 

2 

20060616 

1 

20060725 

1 

20060726 

2 

20060727 

2 

20060728 

1 

20060911 

1 

20060912 

1 

1846Z 


Z 


529Z 

804Z 


17Q5Z 

1836Z 

2046Z 


259Z 

532Z 


1431Z 

1842Z 

2002Z 


54QZ 

739Z 

854Z 


1731Z 


1817Z 


539Z 

825Z 


44QZ 

638Z 


434Z 

836Z 


339Z 

450Z 


348Z 

830Z 


342Z 

437Z 


1438Z 


1533Z 


241Z 

432Z 


3Q4Z 

426Z 


1457Z 


1444Z 


1256Z 
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20060914 


2 


20060915  2 

20061016  3 

20061017  3 

20061018  2 

20061019  3 

20061020  2 

20061127  4 

20061128  2 

20061129  _ 1_ 

20061130  3 

20061201  3 

20061204  _ 1_ 

20061205  3 
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